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Abstract 

These lectures present the analysis of stability and control of long time 
behavior of PDE models described by nonlinear evolutions of hyperbolic 
type. Specific examples of the models under consideration include: (i) 
nonlinear systems of dynamic elasticity: von Karman systems, Berger's 
equations, Kirchhoff - Boussinesq equations, nonlinear waves (ii) nonlin- 
ear flow - structure and fluid - structure interactions, (iii) and nonlinear 
thermo-elasticity. A characteristic feature of the models under considera- 
tion is criticality or super-criticality of sources (with respect to Sobolev's 
embeddings) along with super-criticality of damping mechanisms which, 
in addition, may be also geometrically constrained. 

Our aim is to present several methods relying on cancelations, har- 
monic analysis and geometric analysis, which enable to handle criticality 
and also super-criticality in both sources and the damping of the under- 
lined nonlinear PDE. It turns out that if carefully analyzed the nonlin- 
earity can be taken "advantage of in order to produce implementable 
damping mechanism. 

Another goal of these lectures is the understanding of control mecha- 
nisms which are geometrically constrained. The final task boils down to 
showing that appropriately damped system is "quasi-stable" in the sense 
that any two trajectories approach each other exponentially fast up to a 
compact term which can grow in time. Showing this property- formulated 
as quasi-stability estimate -is the key and technically demanding issue 
that requires suitable tools. These include: weighted energy inequalities, 
compensated compactness, Carleman's estimates and some elements of 
microlocal analysis. 
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1 Introduction 

These lectures are devoted to the analysis of stability and control of long time 
behavior of PDF models described by nonlinear evolutions of hyperbolic type. 
Specific examples of the models under consideration include: (i) nonlinear sys- 
tems of dynamic elasticity: von Karman systems, Berger's equations, Kirchhoff 
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- Boussinesq equations, nonlinear waves (ii) nonlinear flow - structure and fluid 

- structure interactions, (iii) and nonlinear thermo-elasticity. A goal to accom- 
plish is to reduce the asymptotic behavior of the dynamics to a tractable finite 
dimensional and possibly smooth sets. This type of results beside having in- 
terest in its own within the realm of dynamical systems, are fundamental for 
control theory where finite dimensional control theory can be used in order to 
forge a desired outcome for the dynamics evolving in the attractor. 

A characteristic feature of the models under consideration is criticality or 
super-criticality of sources (with respect to Sobolev's embeddings) along with 
super-criticality of damping mechanisms which may be also geometrically con- 
strained. This means the actuation takes place on a "small" sub- region only. 
Super-criticality of the damping is often a consequence of the "rough" behavior 
of nonlinear sources in the equation. Controlling supercritical potential energy 
may require a calibrated nonlinear damping that is also supercritical. On the 
other hand super-linearity of the potential energy provides beneficial effect on 
the long time boundcdncss of semigroups. From this point of view, the non- 
linearity does help controlling the system but, at the same time, it also does 
raise a long list of mathematical issues starting with a fundamental question of 
uniqueness and continuous dependence of solutions with respect to the given 
(finite energy) data. It is known that solutions to these problems can not be 
handled by standard nonlinear analysis-PDE techniques. 

The aim of these lectures is to present several methods of nonlinear PDE 
which include cancelations, harmonic analysis and geometric methods which 
enable to handle criticality and also super-criticality in both sources and the 
damping. It turns out that if carefully analyzed the nonlinearity can be taken 
"advantage of in order to produce implementable control algorithms. 

Another aspects that will be considered is the understanding of control mech- 
anisms which are geometrically constrained. Here one would like to use minimal 
sensing and minimal actuating (geometrically) in order to achieve the prescribed 
goal. This is indeed possible, however analytical methods used are more sub- 
tle. The final task boils down to showing that appropriately damped system is 
"quasi-stable" in the sense that any two trajectory approach each other expo- 
nentially fast up to modulo a compact term which can grow in time. Showing 
this property- formulated as quasi-stability estimate -is the key and technically 
demanding issue that requires suitable tools. These include: weighted energy in- 
equalities, compensated compactness, Carleman's estimates and some elements 
of microlocal analysis. 

The lecture are organized as follows. 

• We start in section 2 with description of main PDE models such as wave 
and plate equations with both interior and boundary damping. Instead of 
striving for the most general formulation, we provide canonical models in 
the simplest possible form which however retains the main features of the 
problems studied. 

• Section 3 deals with wellposedness of weak solutions corresponding to 
these models. 
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• Section 4 describes general and abstract tools used for proving existence of 
attractors and also properties of attractors such as structure, dimension- 
ality and smoothness. Here we emphasize methods which allow to deal 
within non-compact environment typical when one deals with hyperbolic 
like dynamics and critical sources. Specifically a simple but very useful 
method of energy relation due to J. Ball (see [Bal04] and also [MRW98]) 
is presented (allowing to deal with supercritical sources), a version of 
"compensated compactness method" introduced first by A. Khanmame- 
dov [Kha06] (allowing to deal with some critical sources) and a method 
based on " quasi-stability" estimate originated in the authors work [CL04a] 
(and further developed in [CL08a]) which gives in one shot several prop- 
erties of the attractor such as smoothness, finite-dimensionality and also 
limiting properties for the family of attractors. 

• Section 5 demonstrates how the abstract methods can be applied to the 
problems of interest. Clearly, the trust of the arguments is in derivation 
of appropriate inequalities. While wc provide the basic insight into the 
arguments, the details of calculations are often referred to the literature. 
This way the interested reader will be able to find a complete justification 
of the claims made. We should also emphasize that presentation of our 
results is focused on main features of the dynamics and not necessarily on a 
full generality. However, subsections such as Generalizations or Extensions 
provide information on possible generalizations with a well documented 
literature and citations. 

• Section 6 is devoted to other models such as: structurally damped plates, 
fluid-structures, fluid-flow interactions, thermoelastic interactions, Midlin- 
Timoshenko beams and plates. Quantum Zakharov system and Shro- 
dinger-Boussinesq equations. These equations exemplify models where 
general abstract tools for the treatment of long behavior, presented in 
section 4 apply. Due to space limitations, the analysis here is brief with 
details deferred to the literature. 

• Each section concludes with examples of further generalizations-extensions 
and also with a list of open problems. 

Basic notations: Let ficK'^, (i = 2,3, bca bounded domain with a 
smooth boundary F = dft. We denote by H'^{il) i2-bascd Sobolcv space of 
the order s endowed with the norm ||u||s,n = ||''i||Hs(o) a-nd the scalar product 
{u,v)s.n- As usual for the closure of C^(f2) in _ff''(fi) we use the notation 
Hq^H). Below we use the notations: 



We also denote by C{0,T,Y) a space of strongly continuous functions on the 
interval [0, T] with values in a Banach space Y. In the case when we deal with 
weakly continuous functions we use the notation Ct„(0,T, y). 





uvdT, 



'T = 



Ox(o,r), ST = rx(o,T). 
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2 Description of the PDE models 



Wc consider second order PDE models of evolutions with nonlinear (velocity 
type) feedback control represented by monotone, continuous functions, and 
critical-supercritical sources defined on $7 C M.'^, d = 2,3. In what follows, 
we describe models of wave and plate equations with a dissipation occurring 
either in the interior of the domain or on the boundary. Three examples given 
below constitute canonical models which admit vast array of generalizations, 
some of which are discussed in Section 6. We begin with the interior damping 
first. 

2.1 Waves and plates with nonlinear interior damping 
and critical-supercritical sources 

(A) Wave equation with nonlinear damping-source interaction. In a 

bounded domain f2 C K"^ we consider the following wave equation with the 
Dirichlet boundary conditions: 

wtt- Aw + a{x)g{wt) ^ f{w) in ft x {0,T); w = on rx(0,r), (2.1) 

where T > may be finite or oo. We suppose that the damping has the 
structure g{s) = gis+ \s\"^''~^s for some m > 1 and assume the following typical 
behavior for the source f{w) ^ —\w\^~^w, where either 1 < p < 3 or else 

3 < p < min|5, ;^^|. The nonnegative function a € C^(fi) represents the 
support and intensity of the damping. The associated energy function has the 
form: 

E{t) ^ \\\wm? + ^l|V^(t)||^ - f{w{t))dx, (2.2) 

where / denotes the antiderivative of /, i.e. /' = /. The energy balance relation 
for this model has the form 

E(t) + ^ {{ag(wt),wt)))dT = E{s). (2.3) 

When p < 3 the source in the wave equation is up to critical. This is due to 
Sobolev's embedding H^iyi) C L(,{^1) so that f{w) G L2{^) for w of finite en- 
ergy. In this case local (in time) well-posedness holds with no restrictions for m. 
Globality of solutions is guaranteed by suitable a priori bounds resulting from 
either a structure of the source f{w) or from the interplay with the damping, 
see Theorem 3.3 below. 

Long time behavior will be analyzed under the condition of weak degeneracy 
of the damping coefficient a{x) in the so called "double critical case", i.e. p < 3, 
m < 5. We use the compensated compactness criterion in Theorem 4.3 for 
proving existence of an attractor and quasi-stability inequality (4.9) required by 
Definition 4.13 in order to establish both smoothness and finite dimensionality 
of global attractor. 
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When 3 < p < 5 the damping required for uniqueness of solutions needs 
to be superlinear. With sufhciently superlinear damping we show that the cor- 
responding dynamical system is well posed on the finite energy phase space. 
However, the issue of long time behavior, in that case, is still an open problem. 
The known frameworks for studying supercritical sources fail due to nonlinear- 
ity of the damping. More details will be given later. The ranges oi p G [5,6) 
arc considered in [BL08a]. 

We also refer to [CV02, Chapter XIV] where the case of a linear damping 
and a supercrtical force is considered from point of view of trajectory attractors. 

(B) Von Karman plate equation with nonlinear interior damping. Let 

C and a € [0, 1]. We denote Ma = I — a A and consider the equation 

MaWtt + A'^w + a{x)[g{wf) - adivGiVwt)] = [J'{w),w] + P{w) (2.4) 
in f2 X (0, oo) with the clamped boundary conditions: 

d 

w=—w:^0 on rx(0,oo), (2.5) 
on 

where the Airy stress function T{w) solves the elliptic problem 

d 

A^F{w) = -[w.wVinn, with J" = — J" = on T, (2.6) 

on 

and the von Karman bracket [u, v\ is given by 

[u, v] = dl^u ■ dl^v + dl^u ■ dl^v - 2 • dl^^^u ■ d^^^^v, (2.7) 
The damping functions 5 : M M.+ and G : have the following form 

gis) = gis+\sr-^s and G(s, a) = Gi • (s; a) + (Isp-is; a™^ V) (2.8) 

where gi and Gi are nonnegative constants. 

The source term P is assumed locally Lipschitz operator acting from H^{Q,) 
into L2{^) when a = and from H^{il) into H^^{n) when a > 0. 

The associated energy function and the energy balance relation have the form 

Sit) EE i (ww^m' + a\\vw,m') + IwAwitw + IwA^y^w m 

and ^ ^ 

£it)+ [ D{wt)dT = E{s) + ( {{P[w),wt))dT, (2.10) 

J s J s 

where the interior damping form D{wt) is given by 

D{wt) = iiagiwt),wt)) + a{iaG{Vwt),Vwt)). (2.11) 

The case when a > is subcritical with respect to the Airy source. The 
well-posedness and existence of attractors in this case is more standard [CLIO]. 
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When a = the Airy stress source is critical. Wellposedness will be achieved by 
displaying hidden regularity (see Lemma 3.9 below) for Airy's stress function. 
Existence of attractors along with their smoothness will be shown by means of 
quasi- stability inequality (sec (4.9) in Definition 4.13) which can be proved as 
long as P is subcritical. 

(C) KirchhofT-Boussinesq plate with interior damping. With notations 
as in the case of the von Karman plate we consider 

MaWtt + A'^iu + a{x) [g{wt) - a div G(Vwt)] = div [iVw^Vw] + P{w) (2.12) 

with the clamped boundary conditions (2.5). The damping functions g : M i— > 
R+, G : i-^ and the source P are the same as above. The associated 
energy function has the form 

£{t) EE i ilMW + aWVwdtW) + l\\Awit)\\^ + \Vw{x,ttdx (2.13) 

The corresponding energy balance relation is given by (2.10). 

The well-posedness for the case a > is standard. This is due to the fact 
that divlVwpVw € if~^(ri) for finite energy solutions w. The case a = is 
subtle. Its analysis requires special consideration and depends on linearity of 
the damping. Similarly, long time behavior dealt with by using Ball's method 
presented in an abstract version in Theorem 4.4 also requires linear damping. 

2.2 Nonlinear waves and plates with geometrically 

constrained damping and critical-supercritical sources 

This class contains models in the situation when the support of the damping 
a{x) is strictly contained in £7. This is to say suppa(x) C f2o C f2. In addition, 
we considered singular case of interior localized damping which is the boundary 
damping. These models are described next. 

(A) Wave equation with nonlinear boundary damping-source interac- 
tion. In a bounded domain 51 C M'^ we consider the following wave equation 

wtt - Au. + a{x)g{wt) = f{w) in ilx (0, T), (2.14) 

with the Neumann boundary conditions: 

—w + go{wt)^h{w) in rx(0,r). (2.15) 
on 

Concerning the internal damping g and source / terms our hypotheses are the 
same as in Section 2.1(A). The internal damping coefficient a{x) satisfies rela- 
tions 

a e C^{il), a{x) > in 17 with suppa C fio CC 51 

The boundary damping has the form go{s) = g2S + \s\''~^s with q > 1 and 
ff2 > 0. The boundary source has the behavior h{w) ^ ~\w\''~^w, where 1 < 

k < max|3, 
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The associated energy function is given by 

E{t)^]^\\wt{t)f + \\\vwm^^ jji^im^- j^k^m^ (2.16) 

where as above / and h denote the antiderivative of / and h. The energy balance 
relation in this case has the form 

£{t)+ f [iiag{wt),wt))+«goiwt),wt »]dT = £{s). (2.17) 

J s 

Local and global well-posedness result will be presented for the boundary- 
source and damping model. One could consider larger range of boundary and 
interior sources k G [1,4), p € [1,6) [BL08a]. In this case, however, potential 
energy corresponding to the sources is not well defined. Well-posedness result re- 
quires higher intcgrability of initial data [BL08a, BLIO]. More general structure 
of sources and damping can be considered. However, the required polynomial 
bounds are critical. 

In the triple critical case {p = 3,q = 3,ni = 5} with k ~ 2 the existence of 
a global attractor is shown by taking advantage of compensated compactness 
result in Theorem 4.3 (see [CL07a]). With the additional hypotheses imposed 
on the damping smoothness finite-dimensionality of attractor is established by 
proving that the system is quasi-stable (in the sense of Definition 4.13). 

In the case when only boundary damping is active (or internal damping 
is localized) existence and smoothness of attractor require additional growth 
condition restrictions imposed on the damping. This is to say, we need to 
assume p < 3 and 9 < 1 [CLT08, CLT09]. 

(B) Von Karman plate equation with nonlinear boundary damping. 

In a bounded domain fi C we consider the equation 

MaWtt + A'^w + a{x)[g{wt) -adiYGiVwt)] ^ [J'{w),w] + P{w) (2.18) 
with the hinged dissipative boimdary conditions: 

w = 0, Aw = -goi-^wt), onrx(0,r) (2.19) 
on 

Here we use the same notations as in Section 2.1(B). In particular, Airy's stress 
function J^{w) solves (2.6) and the internal damping functions g and G have 
the form (2.8). Concerning the boundary damping we assume that ga^s) ~ 
52 s + |s|'^^s, q > 1. The associated energy function £ has the same form as in 
Section 2.1(B), see (2.9). The corresponding energy balance relation reads as 
follows 



£it) 



dwt \ dw; 



dn dn 

where the internal damping term D{wt) is given by (2.11 



dT = £is)+ I {{P{w),Wt))dT, 

(2.20) 
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Below we show the existence of a continuous semiflow corresponding to (2.18) 
and (2.19). The existence and smoothness of an attractor with boundary damp- 
ing alone is shown under additional hypotheses restricting the growth of the 
boundary damping. In the case of boundary damping one could consider damp- 
ing affecting other boundary conditions such as free and simply supported, see 
[CLIO]. 

(C) Kirchhoff -Boussinesq plate with boundary damping. With the 
same notations as above in a domain C K.^ we consider the equation 

MaWtt + ^"^w + a{x)[g{wt) - a div GiVwt)] = div[|Vu;pVw] + P{w) (2.21) 

with the hinged dissipative boundary conditions (2.19). The internal damping 
functions g and G satisfy (2.8). The boundary damping is the same as the 
previous case, i.e., gois) ^ g2S + |s|''~^s, g > 1. The source P is same as above. 
The associated energy function has the form (2.13). The energy balance relation 
is given by (2.20). 

When a > the wellposedness of finite energy solutions with internal damp- 
ing follows from the observation that the model can be represented as a locally 
Lipschitz perturbation of monotone operator. This approach is also adaptable 
to boundary damping, like in the case of von Karman plate. In the case a = 
the problem is much more delicate. While existence of finite energy solutions can 
be established with general form of the interior-boundary damping, the unique- 
ness of solutions requires the linearity of the damping. In addition, continuous 
dependence on the data depends on time reversibility of dynamics which, in 
turn, does not allow for boundary damping. In view of the above most of the 
questions asked are open in the presence of boundary damping. 

3 Well-posedness and generation of continuous 
flows 

In this section we present several methods which enable to deal with Hadamard 
well-posedness of PDE equations with supercritical Sobolev exponents. By su- 
percritical, we mean sources that may not belong to finite kinetic energy space 
for trajectories of finite energy. While existence of finite energy solutions is 
often handled by various variants of Galerkin method, it is the uniqueness and 
continuous dependence on the data that is problematic. There is no unified 
theory for the treatment of such problems, however there exist methods that 
are applicable to classes of these problems. 

In this section we concentrate on three models described in Section 2.1, where 
each of the model displays difi^erent characteristics and associated difficulties. 
In order to cope with this, different methods need to be applied. We provide 
a show case for the following methods which are critical for proving Hadamard 
well-posedness in the examples cited: 

• Interaction of superlinear sources via the damping. 
Illustration: wave equation. 
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• Cancelations-harmonic analysis and microlocal analysis methods. 
Illustration: Von Karman equations. 

• Sharp control of Sobolcv's embeddings and duality scaling. 
Illustration: Kirchoff-Boussinesq plate. 

1. In the first example the equation can be viewed as a perturbation of a 
monotone operator. However, the superlinearity of the source destroys locally 
Lipschitz character of the semilinear equation. In order to offset the difficulty, 
superlinear damping g(wt) is used. The interplay between the superlinearity 
of potential and kinetic energy lies in the heart of the problem. The presence 
of this damping becomes critical in establishing well defined dynamical system 
evolving on a finite energy phase space. 

2. In the second example, classical regularity results for von Karman nonlin- 
earity, when a = 0, fail to show that the perturbation of monotone operator 
is locally Lipschitz. However, in this case, more subtle methods based on har- 
monic analysis and compensated compactness allow to show that the Airy's 
stress function has hidden regularity property. This estimate allows to prove 
that, contrary to the original prediction, the semilinear term is locally Lipschitz. 
This allows to prove, again, that the resulting dynamical system is well posed 
on a finite energy space. 

3. In the third example, the source of the difficulties is restoring force 

div [iVwpVw] 

which fails to be locally Lipschitz with respect to finite energy solutions in 
the case a = 0. By using method relying on logarithmic control of Sobolev's 
embeddings and topological shift of the energy ("Sedenko's method") we are 
able to prove that the resulting system with linear damping is well posed on 
finite energy space. Critical role in the argument is played by time reversibility of 
the dynamics and linearity of the damping. This excludes superlinear damping 
and boundary dissipation in the case a = 0. 

Conclusion: These three canonical examples present three different methods 
on how to deal with the loss of local Lipschitz property and still be able to 
obtain well-posed flows defined on finite energy phase space. The methods 
presented are transcendental and applicable to other dynamics displaying similar 
properties. Some of the examples are given in Section 6. 

3.1 Wave equation with a nonlinear interior damping - 
model in (2.1) 

The statements of the results 

With reference to the model (2.1) with a{x) > oq > m the following 
assumption is assumed throughout. 
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Assumption 3.1 1. A scalar function g[s) is assumed to he of the form 
g{s) = gis+g2{s), where gi e M and g2{s) is continuous and nondecreasing 
on R with 52(0) = 

2. The source f{s) is either represented by a function and such that 
\f'{s)\ < C{l + \s\P-^) with 1 <p < 3, or else we have that f e C^(R) and 
(i) there exist mg^,Mg^ > such that the damping function g2 satisfies 
the inequality 

mgjsr+' < g2{s)s < Mjsr+\ \s\ > 1. (3.1) 

for some m > 1, (ii) |/"(s)| < C(1 + |s|p~'^) with p satisfying the following 
compatibility growth condition 

3 < p < 6 and p< (3.2) 

m + 1 

While in the case when the source exponent p < 3 we adopt classical semigroup 
definition of the generalized solution (see, e.g., [CL08a]), for supercritical cases 
p > 3 we provide the following variational definition. 

Definition 3.2 (Weak solution) Let (3.1) and (3.2) be in force. By a weak 
solution 0/(2.1), defined on some interval (0,r) with initial data {wo;wi) from 
Hi{fl) X £2(11), we mean a function u e Cw{OjT; H^(fl)), such that 

1. ut e L„+i((0,T) X n)r\C^{Q,T-L2{n))- 

2. For all(j}eC{0,T,H^{n))nC^0,T;L2{n))nLrn+i{[0,T]xn), we have 
that 

{-ut(t>t + VuV(?i) dfldt + [ [ g{ut)(t) dfldt 

Jo Jn 



ut(pdfl 



T 







+ / / f{u)(t) dVldt. (3.3) 



3. lim(ii(t) - u^,(j})isi = and \\m{{ut{t) ~ ui,0)) = for all (j) e i^o(^^). 

Theorem 3.3 Let Assumption 3.1 be in force. We assume that initial data 
{wq^Wi) possess the properties 

w{0) = Wo e -ffo(rj), wt{0) =wie L2{Vl) 

and for p > 5, Wq & Lr{i^) with r = |(p — 1). 
Then 

• There exist a unique, local (in time) solution w{t) of finite energy. This 
is to say: there exists T > such that 

we C([0,T];i/i(r!)), wt e C([0,T];i2(f^)), 

where w is a generalized solution when p < 3 and w is a weak solution 
when condition (3.1) holds with some m > 1 (which is the case when 

P>V- 
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• In this case 3 < p < 5 weak solutions satisfy the energy relation in (2.3). 

• If P 1^ (3.1) holds with some m, then generalized solution satisfies 
also the variational form (3.3) with the test functions (j) from the class 

C(0, T; H^in)) n C\0, T; L2{n)) n Loo(Qt). 

Moreover, the said solutions are continuously dependent on the initial 
data. 

• When p < m the obtained solutions are global, i.e., T = oo. The same 
holds under dissipativity condition: 

- f(s] 

liminf^^>-Ai (3.4) 

>oo S 

where Ai is the first eigenvalue of —A with the zero Dirichlet data. 

• When p > 3, p > m and f{s) = |s|^^^s local solution has finite time 
blow-up for negative energy initial data. 

Remark 3.1 In the case when p £ [1, 3] one constructs solutions as the Umits 
of strong semigroup solutions [CL07a]. These are referred to as generalized 
solutions[Ch07a,, CL08a] . It can be shown that generalized semigroup solutions 
do satisfy the variational form in question and variational weak solutions are 
unique [CL08a]. In order to justify this, it suffices to consider only the nonlinear 
damping term 17(5) with gi — 0, since the latter provides a bounded linear 
perturbation-thus it does not affect the limit passage. The energy inequality 
provides Li bound for g{wn). This allows to apply Dunford-Pettis compactness 
criterion (see later) in order to transit with the weak limit in Li, and thus 
obtaining variational form of solutions (3.3) with appropriate test function as 
specified in the third part of Theorem 3.3. 

Sketch of the proof 

Without loss of generality wc may assume a{x) = 1 and 171 = 0, so 32(5) ~ gis) 
(these parameters have no impact on the proof). 

Case 1: Critical/subcritical sources and arbitrary monotone damping. 

In this case the proof is standard and follows from monotone operator theory. 
Setting the second order equation as a first order system 

Wt+A(W)=F(W), t>0, W{0) ^ {wo;wi) (3.5) 

on the space H = i?Q (51) x ^2(57) where W = {w; wt) and 

allows to the use of monotone operator theory (see, e.g., [Bar 76, Sho97] for the 
basic theory). Indeed, A is maximally monotone and F{W) is locally Lipschitz 
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on H. This last statement is due to the restriction p < 3 and Sobolev's embed- 
ding ffi(r^) C LQ{n). General theorem (see [CEL02] and also [CLIO, Chapter 
2]) allows to conclude local unique existence of semigroup solutions. 

The arguments leading to variational characterization of weak solutions are 
given in the Appendix [CL07a] . It can be based on weak Li compactness criteria 
due to Danford and Pettis (see [DS58, Section IV. 8]) which states that the set 
M C Li{Q) is relatively compact with respect to weak topology if and only if 
this set is bounded and uniformly absolutely continuous, i.e.. 



Ve > 3(5 > : mes{E) < S 



sup 

feM 



fdQ 



< e. 



We want to show that M = {g(?ii„)} is weakly Li{Q) compact for a sequence 
of strong solutions which also converges almost everywhere. For this using 
dissipation inequality g{wn)wndQ < C which (by splitting the integration 
into \wn\ < R and > R) implies that 

/ \g{wn)\dQ < 4 / 9{wn)wndQ + ffflmes {E), \/R>0. 
JE R Jq 

Case 2: supercritical source. In this case superlinearity of the damping 
is exploited. The key role is played by energy inequality obtained first for 
smooth solutions corresponding to truncated problem where / is approximated 
by locally Lipschitz function fx defined in [BLIO]. Locally Lipschitz theory 
(Case 1 above) allows to obtain the energy estimate 



{{g{wt),Wt))ds<E{Q)+ (ifK{w),wt))ds (3.6) 



with E{t) = i(||u;f(i)|p + ||Vw(t)|p). Exploiting the growth conditions imposed 
on g and / yields 



E{t) + \\wt 



im+l 



< C 



£;(0) + |1/k(u;)|1- 



and Sobolev's embeddings iJ^(il) C L^(^) along with the condition 
p(m + 1) < 6™ implies 



E{t) 



E{Q)^L\ t 



The above estimate followed by a rather technical limit argument (with K — > oo) 
[BLIO] allows to establish local existence of weak solutions. 

Energy identity. This is an important step of the argument. By using finite 
difference approximation for the velocity Wtit) one shows that weak solutions 
satisfy not only energy inequality but also energy identity. The detailed argu- 
ment is given in the proof of Lemma 2.3 [BLOSa]. 
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Uniqueness of the solutions is more subtle. We need to show that any two 
solutions from the existence class, i.e. satisfying 



is uniquely determined by the initial data. To this end let's consider the differ- 
ence of two solutions: z = w — u, where both w and u are finite energy weak 
solutions specified as above with the same initial data (wq] wi). Denoting by 



owing to energy inequality satisfied for weak solutions and using bounds imposed 
on the source function /(s) one obtains after some calculations: 



The following lemma is critical: 

Lemma 3.4 Letp<5. ThenVe > andy{wo;wi) G Hq^H) x L2{^) such that 
||wo||i,n + llwillo.n < R, there exist a constant Ce(R,T) > such that 



This lemma is specialized to the case when p < 5. For p € [5,6) there is another 
term in the inequality above which accounts for the fact that the potential 
energy term may not be in ii. In order to focus presentation wc omit this term 
and refer the reader to the source [BL08a] . 

Proof. Using integration by parts we have that 



w £ cU[o,T],H^{n)nci{[o,T],L2in)), wt e l,„+i(Qt) 



E,{t) = \\\/z\\^ + \\zt\\^ 




(3.7) 



\Rf{T)\<eE,{T) + C,{R,T) TE,{T) 





{{f{w)^f{u),z)) 



2((/'H,~-')) 



J 
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It is convenient to use the following elementary estimate valid with any finite 
energy element z such that z(0) = 0: 

f ft 2 /■* /■* 

\WW ^ / / zt{s)ds dx<ti \\zt{s)\\'^ds <2t I E,{s)ds. (3.8) 



Estimate for J^^ z"^ {t)\ut{t)\ dxdt. Wc use Holder's Inequality with p = 3 
and (7 = 3/2 and obtain: 

z\t)\ut{t)\ dQ< [ ||z(t)|li^(j,) • |kt(t)||L3/.(f2) dt (3.9) 
Jn Jo 

This leads to 

r [ z\t)[\wt{t)\ + \ut{t)\] dQ < C[R) r E,{t)dt. (3.10) 
Jo Jn Jo 

Estimate for \u{T)\p~^z'^{T) dx. First, we write 

\u{T)\P~^z\T) dx < \\z{T)f + f \uiT)\P-^z\T) dx. (3.11) 

Jf2n{|u{T)|>l} 

The first term is estimated in (3.8). The argument for the second term is given 
below. We consider here the supercritical case 3 < p < 5. The corresponding 
estimate for p £ [5,6) is given in [BL08a]. 

Since \u{T)\ > 1, there exists eo > such that \u{T)\p^'^ < |u(T)|'*-"«. We 
choose e < ^ and apply Holder's inequality with q = and q ~ 2{i-e) ■ 
then use Sobolev's embeddings, the fact that (4 — eo) 2(1-6) — ^' interpolation 
and Young's inequality to obtain 



/ \uiT)\P-^z^{T) dx 

Jnn{\u{T)\>i} 



< 



l + 2e 2(l-g) 



( 1^ |z(r)|iT37 dx) ' ( j \u{T)f-^ dx) ' (3.12) 



< C\\z{T)rH,^,\\u{T)\\['^^^^< eE^{T) + C,{R)\\z{T)\\\ 

2(l-e) 

Combining (3.11) with (3.8) and (3.12), wc obtain the final estimate in this case: 
\u{T)\P-'^z^{T) dx < eE,{T) + C,{R)T [ E,{t) dt (3.13) 



Estimate for J^\u{t)\P~'^\ut{t)\z'^(t) dQ. Since for \u{t)\ < 1, we obtain 
the term estimated in the previous case, it is sufficient to look at 

/ / \u{t)\p-^\utit)\z'^{t) dQ with n^nn {\u{t)\ > i}. 

Jo Jn 
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We start by applying Holder's inequality with g = 3 and q = 3/2: 



We use Holder's inequality again, with q = and q = tA- and notice that 



(3.14) 

^1 I 1 1 I II o n /H n 1^+' T 

p-2 ^ y 6~p ' 



< m + 1 ^ p < Hence (3.14) becomes: 

z.(t)r2|^.,(t)|z2(i)<c ri?.(i)||u(t)||^;(%h,(t)IU , (o) 



Jn Jo 



(3.15) 



<c / i?,(oik(i)ii^;(%ikt(t)iiw(o). 



Since ||M(t)||i,si < Cij_T we obtain: 

\u{t)\P~^\ut{t)\z\t) dQ < Cr^t f E,{t)\\ut{t)U,„^^^n) dt. (3.16) 
Jn Jo 

Similarly, following the strategy used in the previous step, we obtain the estimate 
\u{t)\P-'\wt{t)\z^{t) dQ < Cr,t r E,{t)\\wt{t)\\L^^^(^) dt, 

JO„ JO 

^ / \w{t)\P-^\ut{t)\z\t) dQ < Cr.t r EM\ut{t)\L^^,in) dt 
Jn^, Jo 

with (iw = n {|w(i)| > 1}. This completes the proof of Lemma 3.4. ■ 
Completion of the proof: In (3.7), we use Lemma 3.4 to obtain: 

E, (T) < eE,{T) + C,{R,T)[ TE, (T) 

rT (3.17) 

+ / (1 + ll"t(i)IU,„+i(o) + ||«t(i)IU,„+i(o)) E,{t) dt 
Jo J 

for every e > 0. In (3.17), we choose e and T such that e + C,{R,T)T < 1/2 
(for uniqueness it is enough to look at a small interval for T, since the process 
can be reiterated) and apply Gronwall's inequality to obtain that Ez (T) = for 
all T < Tmax, where [0, Tmax] is a common existence interval for both solutions 
w and u. 

The above argument completes the proof of uniqueness of finite energy solu- 
tions. This result along with energy equality leads to Hadamard wellposedness. 
Details are in [BL08a, BLIO]. 

Finite-time blow up of solutions is established in [BL08b]. This is accomplished 
by showing that appropriately constructed anti-Lyapunov function [STV03, 
Vit99] blows up in a finite time. 
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Generalizations-Extensions 

• In the subcritical and critical case additional regularity of weak solutions 
equipped with more regular initial data can be established. Quantitative 
statements are in [CL07a]. 

• The wave model discussed is equipped with zero Dirichlct data. However, 
the same result holds for Neumann or Robin boundary conditions. 

• In the subcritical and critical case, p < i, the same results hold with 
localized damping. This is to say when suppa(a;) C Oo C O. 

• The analysis of wellposedness for the same model in the case when solu- 
tions are confined to a potential well are given in [BRTll]. 

Open question: Obtain the same result with partially localized damping 
suppa(x) C 51o C in the supercritical case p > 3. 

3.2 Von Karman equation with interior damping - model 
(2.4) 

The statement of the results 

With reference to the model (2.4) the following assumption is assumed through- 
out. 

Assumption 3.5 1. Scalar function g{s) is assumed to be continuous and 
monotone on M with g{0) ~ 0. 

2. In the case a > we assume that the damping G has the form G{s, a) = 
(5i(s); (72(c)), where gi{s), i = 1,2, are continuous and monotone on R 
with gi(0) = 0. Moreover, they are of polynomial growth, i.e., \gi{s)\ < 
C(l + \s\''-'^) for some q>l. 

3. The source P{w) is assumed locally Lipschitz from H^{VL) into [Ha{^)]' , 
where 

U^m = {%f^, « ca.e a = 0; 

[ ^o(") the case a > 0. 

Below we also use the following refinement of the hypothesis concerning the 
source P. 

Assumption 3.6 We assume that P{w) = —Pq{w) + Pi{w), where 

1. Pa{w) is a Frechet derivative of the functional Ilo{w) and the property 
holds: there exist a,b € M. and e > such that 

IloH+a\\w\\l_^^,,>b, {{Po{w),w))+a\\w\\l_^^,,>b, e H^{n); 

2. there exists K > and £ > such that \\Pi{w)\\ < K\\w\\2~e,n for all 
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Remark 3.2 A specific choice of interest in applications is P{w) = [Fo,w] —p, 
where [•, •] denotes the von Karman bracket (see (2.7)), Fq G H^+^ (n) n H^{n), 
^ > 0, and p e 1/2 (f^)- The term Fq models in- plane forces in the plate and p is 
a transversal force. The given source complies with all the hypotheses stated. 
In the case a > the hypotheses concerning P{w), Fq and p can be relaxed. 
However we do not pursue this generalizations and refer to [CLIO]. 

In addition to the concept of generalized (semigroup) solutions [CL08a, CLIO] 
we can also define weak solutions. 

Definition 3.7 (Weak solution) By a weak solution of (2.4), defined on an 
interval [0,T], with initial data (uq]Ui) we mean a function 

ueC^,iO,T;H^in)), ut e a40,T;Ha(r!)) 

such that 

1. For all <f) e H^{n), 

{{ut{t), 0)) + a{{\/ut{t), ^cj))) + [ {{An, A0)) dt (3.19) 

Jo 

+ f [{i9{ut),(l>)) + a{{G{Vut{t))y<j>j)] dt 
Jo 

= ((ui»)+a((Vui,V0))+ / iiPiu) + [Fiu),u],(j>)) dt. 

Jo 

2. ]ha{u{t) - Mo, 0)2,0 = 0. 

Here as above Cn]{0,T,Y) denotes a space of weakly continuous functions with 
values in a Banach space Y . 

Theorem 3.8 Under the assumption 3.5 for all initial data 

wiO) = wo e H^in), wtiO) ^wi eH4n) 

(with Hai^) defined by (3.18)j there exist a unique, local (in time) generalized 
(semigroup) solution of finite energy, i.e., there exists T > such that 

w e C{[0,T];H^{n)),wt G C([0,T];H„(17)). 

Moreover, 

• //, in addition in the case a — the damping g(s) is of some polynomial 
growth, generalized solution becomes weak solution. A weak solution is 
also continuously dependent of the initial data. 

• The solutions are global provided Assumption 3.6 holds. In this case a 
bound for the energy of solution is independent of time horizon. 
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Sketch of the proof 

We concentrate on a more challenging case when a = 0. The case a > can be 
found in [CLIO, Chapter 3]. 

For the proof of this result the support of the damping described by a{x) 
plays no role. Thus, without loss of generality we may assume that a{x) = 1. 

The following two lemmas describing properties of von Karman brackets can 
be found in [CLIO]. The first Lemma is critical for the proof of wcU-posedness 
and the second Lemma implies boundedncss of solutions that is independent on 
time horizon. 

Lemma 3.9 Let A^^ denotes the map defined by z = A~^/ iff 

d 

A^z — f in il and f = —f — on F. 

on 

Then 

\\A^'^[u,v]\\w2,oa(^a) < C\\u\\2A-i\\v\\2Ai, 

where the von Karman bracket [w, u] is defined by (2.7). In particular, for the 
Airy stress function T we have the estimate \\J-'{w)\\w^^°°{n) ^ C'll'^^ll? n- 

Notice that standard result [Lio69] gives 

\\A-^[u,v]h.,,n<C\\u\\2M\M2,n 

which does not imply the result stated in the lemma. 

The second lemma [CLIO] has to do with a control of low frequencies by 
nonlinear term. 

Lemma 3.10 Let u e H'^{Vl) n Hq{VI). Then for every e > there exists 
Me > such that 

Equipped with these two lemma, the proof of global well-posedness follows stan- 
dard by now procedure: 

Step 1: Establish maximal monotonicity of the operator 

A=( ^ ^ 

on the space H = -ffo(^) ^ ^2(^1) with 

D{A) = {{u; v) e H^{n) X L2{n) : A^u + g{v) € L2{^)} 

This follows from a standard argument in monotone operator theory [Bar76, 
Sho97]. 
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step 2: Consider nonlinear term as a perturbation of A: 



F{W) = 





[T{w),w] +P{w) 



Step 3: Show that F{W) is locahy Lipschitz on H. This can be done with the 
help of Lemma 3.9 which imphes the following estimate 



\mu),u] - [Tiw),w]\\ < C(l + \\u\\ln + \\w\\ln)\\u - wh, 



(3.20) 



The above steps lead to well-defined semigroup solutions, defined locally in 
time, which are unique and satisfy local Hadamard well-posedness. 

The final step is global well-posedness for which a priori bounds are handy. 
The needed a priori bound results from the property of von Karman bracket 
described in Lemma 3.10. Indeed, in order to claim global well-posedness, it 
suffices to notice that the relation 



dt 



||A.FH|p + noH 



((PiH,w;t)) 



implies the a priori bounds for the energy function £{t) given by (2.9) with 
a = via energy inequality and Gronwall's lemma. 

Stationary solutions. For description of the structure of the global attractor 
we also need consider stationary solutions to problem (2.4). We define a sta- 
tionary weak solution as a function w G H^{Q.) satisfying in variational sense 
the following equations 



A^w = [J"(w),-u;] + P(u;) in f]. 



d 

— w = on F, 
on 



(3.21) 



where the Airy stress function F{w) solves the elliptic problem (2.6). 



Proposition 3.11 Let Assumption 3.6 be in force. In addition assume that P 
is weakly continuous mapping from Hq{VI) into H~^{fl). Then problem (3.21) 
has a solution. Moreover, the set of all stationary solutions is compact in 

For the proof of this proposition we refer to [CR80], see also [CLIO, CiaOO]. We 
note problem (3.21) may have several solutions [CR80]. However in a generic 
situation the set AC is finite (see a discussion in [CLIO, Chapter 1]). 



Generalizations-Extensions 

1. By assuming more regular initial data one obtains regular solutions. Pre- 
cise quantitative statement of this is given in [CLIO]. Regular solutions 
are global in time. 
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2. Related results hold when a > 0. The analysis here is simpler and there 
is no need for sharp results in Lemma 3.9, see [CLIO]. 

3. One can consider other boundary conditions such as simply supported or 
free or combination thereof, [CLIO]. 

4. Damping can be partially supported in fl. This has no effect on the 
arguments. 

5. A related wellposedness result holds for non-conservative models with ad- 
ditional energy level terms that are non-conservative (for instance Vw ■ ^ 
for some smooth vector ^f). 

Open question 

Full von Karman systems that consist of system of 2D elasticity coupled with 
von Karman equations. This model accounts for in-plane accelerations. As 
such, there is no decomposition using Airy's stress function. The nonlinearities 
entering are super-critical (even in the rotational case when a > 0) and hidden 
regularity of Airy's stress function plays no longer any role. For such model with 
a > one can still prove existence and uniqueness of solutions, by appealing to 
"Sedenko's method" (see [Sed91]). Even more, full Hadamard well-posedness 
and energy identity can also be proved in that case [Las98, KL02]. However, 
the problem is entirely open in the non- rotational case a — 0. In this latter case 
only existence of weak solutions, obtained by Galerkin method, is known. 

3.3 Kirchhoff-Boussinesq model with interior damping - 
model in (2.12) 

As before we concentrate on the most demanding case when a ~ 0. For the case 
a > we refer to [CL08a]. The main challenge of this model is the presence 
of restorative force term div [[Vwl^Vw] which is not in L2{^) for finite energy 
solutions w. This is due to the failure of Sobolev's embedding C Loo in two 
dimensions. We also assume, without loss of generality for the well-posedness, 
that a{x) = 1. 

The statement of the results 

With reference to the model (2.12) with a = the following assumption is 
assumed throughout. 

Assumption 3.12 1. g{s) is continuous and monotone on M. with g{0) = 0. 
In addition g{s) is of polynomial growth at infinity and also a(x) = 1 
(without loss of generality for the well-posedness) . 

2. The source P{w) is assumed locally Lipschitz from H^{Q.) into L2{^). 

Below we deal with weak solutions. 
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Definition 3.13 (Weak solution) By a weak solution of (2.12), with a = 
and a{x) = 1, defined on some interval {0,T), with initial data (mo;ui) we mean 
a function u G (7^,(0, T; ffg (17)) such that Ut S (7^,(0, T; L2(ri)) with g{ut) G 
Li{Qt) and 

1. For all (j) e Hl{n), 

ut{t)(j>dn+ [ [ {AuAcI)) dndt + [ [ g{ut)(j) dQ.dt 
Jo Jn Jo Jn 

= [ ui(j)dn+ [ [ [P(u)0- |Vi«p(Viy, V(/))] dndt (3.22) 
Jn Jo Jn 

2. \im{u{t) ~ u„,(t))2,n = 0. 
t^o 

Theorem 3.14 Let a ~ 0. Under the assumption 3.12 for all initial data 

w{0)=woeH^{n), wt{0) = wi e L2in) 

there exist a local (in time) weak solution of finite energy. This is to say: there 
exists T > such that 

w€C^{[0,T];H^in)), wt€C^i[0,T];L2in)). (3.23) 

Moreover 

• Under additional assumption that g{s) is linear the said solution is unique. 
In addition energy identity is satisfied for all weak solutions which are also 
Hadamard well-posed (locally). 

• The solutions are global under the following dissipativity condition: for 
every 5 > there exists Cs > such that 



iiPiwis)),wtis)))ds < 5£{t) + Cs 







£{0)+ / £{s)ds 







(3.24) 



for any function w(t) possessing the properties in (3.23). 



Remark 3.3 A specific clioice of interest in applications is Boussinesq source 
given by P{w) = A[w^]. This source complies with all the hypotheses stated. 
Indeed, 



dt 



which then gives 



{{Aw\wt)) = -|((|VH^u;)) + ((|VH^^t)). (3.25) 
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Since 



\\w{t)r<\\wiO)r + ^ \\w{T)\\\\wt{T)\\dT<C 







£(0)+ / £{T)dT 







relation (3.25) quickly implies the conclusion desired in (3.24). 
3.4 Sketch of the proof 

As mentioned before the challenge in the proof of Theorem 3.14 is to handle the 
lack of local Lipschitz condition satisfied by the restorative force. 

We explain the main steps. Full details are given in [CLll, CL06b] 

Step 1 (Existence): It follows via Faedo-Galerkin method. This step is 
standard. 

Step 2 (Uniqueness): For this we use Scdenko's method which based on writ- 
ing the difference of two solutions in a split form as projection and coprojection 
on some finite dimensional space and estimating. 
We start with some preliminary facts. 

We introduce the operator A in L2(ri) by the formula Au = A^m with 
the domain D{A) = H'^{il) n Hq{Q). The operator A is a strictly positive 
self-adjoint operator with the compact resolvent. Let {ck} be the orthonormal 
basis in L^(r2) of eigenvectors of the operator A and {X^} be the corresponding 
eigenvalues: 

Ae-k^Xkek, fc = l,2,....; < Ai < A2 < . . . 

W also note that for every s <E [0, 1/2] we have DiA") = Hi^^i^), s ^ 1/8, 3/8. 
Moreover, the corresponding Sobolev norms are equivalent to the graph norms 
of the corresponding fractional powers of A, i.e. 

ci\\A'u\\ < \\u\\4s < C2||^'u||, u e D{A'), (3.26) 

for all admissible s E [0, 1/2]. 

The following assertion is critical for the proof. 

Lemma 3.15 ([BC98]) Let P^v be the projector in L'^in) onto the space span- 
ned by {ei, 62, e^v} and f{x) <E D{A^^'^). Then there exists No > such 
that 

max|(PA,/)(x)| <C-{log(l + Ajv)}^/' || / 111 (3.27) 
for all N > Nq. The constant C does not depend on N . 

Remark 3.4 For the first time a relation similar to (3.27) was used for unique- 
ness in some shell models in [Sed91]. Latter the same method was applied for 
coupled 2D Schrodinger and wave equations [CS05, CSll], for the inertial 2D 
Cahn-Hilliard equation [GSZ09] and also for some models of fluid-shell interac- 
tion [CR12]. One can show (see [CLIO, Appendix A] and the discussion therein) 
that (3.27) is equivalent to some Brcsis-Gallouet [BG80] type inequality. 
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Assume that a = 0, a{x) = 1 and g{s) = ks in (2.12). Let wi{t) and W2(i) 
be two weak solutions of the original problem (2.12) with the same initial data 
and w(t) = wi{t) — W2{t). Then WN{t) = PN'w{t) is a solution of the linear, but 
nonhomogenous problem 

Wtt + kwt+A'^w^{PNM){t), xen,t>0, (3.28) 

with the boundary and initial conditions 

dw 



^ 0, w\t=o = 0, dtw\t=o ^ 0. 
on 



Here 

M{t) = div [\Vwi{t)fVwi{t) - \Vw2{t)\''Vw2it)] + [P{wi{t)) - P{w2{m . 
Multiplying the equation (3.28) by A~^^^Wt and integrating we obtain that 

WA-'^^PNWtim' + \\A'/^PNw{t)W' <cf \\A-^'Hl{T)\\\\A-^'^wt{T)\\d7 

Jo 

for all t G [0,T]. From here after accounting for (3.26) we obtain 

II m{t) \Wn + II Ht) \\1<C- f II M{t) W-is, ■ II Wt{T) dr. 



In particular this implies that 

II w{t) \\i,n<C- / II M{t) W-is-idT. (3.29) 



The inequality above is the basis for further estimates. Below we also use the 
estimate which follows from the definition of weak solutions: 

sup {||wi(t)||2,o + ||w2(t)||2,o}<i?, (3.30) 
te[o,T] 

where i? > is a constant. Using Lemma 3.15 we can estimate the quantity 
||i\/(t)||_i in the following way: 

||A/(t)||-i < Ci • log(l + \n)- II w{t) 111 +C2 ■ X^+t (3-31) 

for some < s < 1 and with the constants Ci and C2 depending on R from 

(3.30) (for details we refer to [CL06b]). Therefore it follows from (3.29) and 

(3.31) that ip{t) = ||ti;(i)||i satisfies the inequality 

m < Ci ■ log(l + \n) / ^{T)dT + C2-T- A^^/f , t e [0, T], 



for some < s < 1. Thus using Gronwall's lemma we conclude that 

m < C2 • r • A^^/f • (1 + XNf'\ t e [o,T]. 
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If we let N ^ oo, then for < t < to = s ■ (4Ci)~^ we obtain '0 = 0. Thus 
wi{t) = W2{t) for < t < to. Now we can reiterate the procedure in order to 
conclude that wi{t) = 'W2{t) for all < i < T", where T is the time of existence. 
This completes the proof of uniqueness. 

Step 3 (Energy identity) : Energy inequality relies on a standard weak lower- 
semicontinuity argument. Since the system is time reversible, one obtains energy 
inequality for the backward problem. Combining the two inequalities: forward 
and backward, leads to energy identity (2.10) with S given by (2.13) satisfied 
by weak solutions. Details arc given in [CLll]. 

Step 4: Equipped with uniqueness and energy identity standard argument 
furnishes Hadamard well-posedness. 

Step 5: The last step is to extend local (in time) solutions to the global ones. 
This is based on a priori bounds postulated by nonlinearitics in (3.24). 

Generalizations- Extensions 

1. Higher regularity of solutions can be proved by assuming more regular 
initial data. Quantitative statements are given in [CLll]. 

2. Rotational models, when a > 0, can be considered without extra difficulty. 
In fact, in this subcritical case one obtains full Hadamard wellposedness 
also in a presence of nonlinear damping subject to the same assumptions 
as in the von Karman case. Some details can be found in [CL08a]. 

3. We can also consider different boundary conditions such as hinged and free. 
Free boundary conditions are most challenging due to intrinsic nonlinearity 
on the boundary, [CL06b]. 

4. More general structures of restoring forces can be also considered, see 
[CL06a, CL06b, CL08a, CLll]. 

5. The support of the damping may be localized to a small (or even empty) 
subset of Vlo C Vl. This will not affect a finite time behavior of solutions. 

Open Problem: Uniqueness of weak solutions with nonlinear damping. 

3.5 Wave equation with boundary source and damping - 
model in (2.14) 

With reference to the model (2.14) and (2.15), where, for simplicity, we take 
g = and / = 0, the following assumption is assumed throughout. 

Assumption 3.16 1. Scalar function go (s) is assumed to be continuous and 
monotone on M with (?o(0) = 0. 
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2. The source h{s) is represented by a function such that \h"{s)\ < C(l + 
Isj*^"^), where 2 < k < 4, and the following growth condition is imposed 
on the damping go{s) with the constants mgg,Mgg > 0; 

mgM'^' < 9ois)s < Mg„\s\'^\ \s\ > 1, with q > ^ > 1. (3.32) 

When the damping is suhlinear i.e., q £ (0,1), then the the source is 
required to satisfy \h'{s)\ < C(l + jsj'^^"'^), 1 < fc < 4 and the condition in 
(3.32) should he satisfied for this q € (0, 1). 

Definition 3.17 (Weak solution) By weak solution of problem (2.14) and 
(2.15) with .9 = and / = 0, defined on some interval (0,T) with initial data 
{uq]Ui), we mean a function u € Cto(0, T; iJ^(r2)) such that 

1. ut e C.^(0,T;L2(r2)) and m e Lg+i{T.T), where St = [0,T] x T, 

2. For all (j) e C(0, T, i?i(17)) n C^{Q,T]L2{n)) f] Lq+i{Y,T), 

{-ut(f)t +VuS/(f>) dVldt+ I / gQ{ut)(f> dndt 



Jn 



ut(j)dfl 



n 



T 







T 







h{u)(l) dndt. (3.33) 



3. lim{u{t) - UQ,(j))i^n = and lim((M((t) - ui,0)) = for all 4> e H'^{n). 

Theorem 3.18 Let f = 0, g = and Assumption 3.16 be in force. Let initial 
data {wq;wi) be such that 

w{0)=woeH\n), wt{0)^wieL2{n), 

and also wo|r € L2fc_2(r) when fc > 3. Then there exist a unique, local (in 
time) weak solution of finite energy. This is to say: there exists T > such that 

w e C{[0,T];H\n)), wt € C{[0,T];L2m). 

Moreover, 

• When k < 3 the energy identity holds for weak solutions and weak solution 
is continuously dependent on the initial data. 

• When k < q the obtained solutions are global, i.e., T = oo. The same 
holds under dissipativity condition: —h{s)s > 0. 

• When l<k<3, k>q and h{s) = jsj'^^^s local solution blows up in a 
finite time for negative energy initial data. 
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Remark 3.5 In the case when h[s) = as, no growth conditions imposed on go 
are required. In fact, in that case the obtained solution is semigroup generalized 
solution. 

When \h"{s)\ < c variational form of the solution can be obtained with more 
relaxed hypotheses imposed on high frequencies of the damping. For instance, 
we can assume 

liminf > 0, |5o(s)| < c[l + \s\^] (3.34) 

s|— fCxD S 

Under the above condition one can show that the generalized solution satis- 
fies also variational form (3.33) with the test functions (f> G C{0,T; H^{il)) H 
C^{0,T;L2{n))nC(T,T) [CLOla] 

Reference to the proofs 

Notice that the damping g^ is assumed active for aU values of the parame- 
ter fc > 1. This is unlike interior source where only supercritical values of 
the source require presence of the damping. In the boundary case, however, 
Lopatinski condition is not satisfied for the Neumann problem and this neces- 
sitates the presence of the damping which, in some sense, forces Lopatinski 
condition [Sak70]. This is manifested by the fact that H^{Q) solutions of wave 
equation with Neumann boundary conditions do not possess H^{12) boundary 
regularity (unlike Dirichlet solutions) . It is the presence of the damping which, 
in some sense, recovers certain amount of boundary regularity [LTOO, Vit02a]. 

The proof of Theorem 3.18 follows similar conceptual lines as in the interior 
case. There arc however few subtle differences due to the presence of undefined 
traces in the equation. In the case when |ft."(s)| < C and ms^ < go{s)s < AIs^ 
the proof is given in [CL07a]. For the remaining values of the parameters 
Hadamard wellposedness is proved in [BLOSa, BLIO]. The analysis of the bound- 
ary case is more demanding than in the interior case. The corresponding differ- 
ences are clearly exposed in the structure of a counterpart to Lemma 3.4, which 
in the boundary case has a different form given in Lemma 4.2 in [BLOSa]. 

The last statement in Theorem 3.18 regarding finite time blow up of energy 
is proved in [BL08b] . 

Generalizations-Extensions 

1. Additional regularity of solutions corresponding to more regular initial 
conditions when q < 3 and |/i"(s)| < C is given in [CL07a]. 

2. One could obtain a version of Theorem 3.18 that incorporates both damp- 
ings: internal and boundary. This is done in [BLOSa, BLIO]. 

3. Well-posedness theory for small data taken from potential well (fc < 3) is 
also available [BRTll, Vit02b]. 

Open questions: 
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1. Well-posedness of finite energy solutions theory without the boundary 
damping when k < 2 (so that h{u) € i2(r) for finite energy solutions. 

2. Interplay between boundary and interior damping. Is it possible to show 
well-posedness of finite energy solutions for supercritical case p > 3 with 
boundary damping only (i.e. g = 0)? 

3. The same question asked for finite time blow results. Does boundary 
source alone leads to blow up of energy in the presence of a boundary 
damping when q < k? 

4. Boundedness of solutions when time goes to infinity and sources are gener- 
ating energy. Quantitative description of the behavior of global solutions 
- when p < m and k < q. 



3.6 Von Karman equation with boundary damping - 
model (2.18) 

With reference to the model (2.18) and (2.19) the following assumption is as- 
sumed throughout. 

Assumption 3.19 1. Scalar function g{s) and gQ{s) are assumed to be con- 
tinuous and monotone on M with at 0. The rotational damping G (the 
case a > 0) satisfies Assumption 3.5(2). 

2. The source P{w) is assumed locally Lipschitz from H^{fl) into H^^(fl) 
when a > and from H^(fl) into L2{i^) when a = (see Remark 3.2 
concerning a specific choice of interest in applications). 

We concentrate on a more challenging case when a = 0. In addition to the 
concept of generalized (semigroup) solutions we also define weak solutions. 

Definition 3.20 (Weak solution) By a weak solution of (2.18) and (2.19) 
with a = and with initial data (uojwi), defined on some interval (0,T), we 
mean a function u G C^{Q,T; H'^{n) r\ Hl{n)) such that ut € C^{Q,T; L2{^)) 
with g(ut) e Li{Qt), go{{d/dn)ut) E Li(I]t) and 

1. For all (f) e ij2(0) n Hl{VL), 

(AuA6)dndT+[ [ g{ut)<b dVldr + I [ go{-^Ut)-^(l)dTdT 
Jn Jo Jn Jo Jr on dn 

{ut{t) - ui)<t)dVl +11 {P{u) + [F{u),u])<t) dVLdr. (3.35) 
Jo Jn 

2. lirn(u(t) — Mq, 0)2,0 = and lini((uf(t) — ui,(/))) = 0. 
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Theorem 3.21 Let a ~ 0. Under Assumption 3.19 for all initial data 

w{Q) = woe H^{n) n H^{n), wt{0) = wi e L2{n) 

there exist a unique, local ( in time ) generalized ( semigroup ) solution w of finite 
energy, i.e., 

3T>0 : we C([0, T]; H\n)), wt G C([0, T]-L2{^)). 

Moreover 

• //, in addition, the damping g{s) and go{s) are of some polynomial growth 
and {go{x) — goiy)){x — y) > c\x — for some r > 1, then a generalized 
solution becomes also weak solution. Weak solution is also continuously 
dependent on the initial data. 

• The solutions are global under Assumption 3.6 concerning P{w) with the 
energy bound independent of the time horizon. 

Reference to the proof 

Well-poscdncss of Von Karman plates with the boundary damping oceurring in 
the moments have been proved in [CLIO]. See also [JL99, HL95]. 

3.7 Generalizations 

1. Additional regularity of finite-time solutions for more regular initial data, 
[CLIO, HL95]. 

2. Similar problems can be formulated for plate equations with the damping 
acting in shears and torques [Lag89] . Thus the boundary conditions under 
consideration are "free" with the feedback control given by go{wt) acting 
on the highest third order boundary conditions [CL04b, CL07b]. 

3. Combination of interior and boundary damping can also be considered by 
combining the methods. 

4. Model with rotational inertial forces a > subject to boundary damping, 
[CLIO, HL95]. 

5. Regular solutions for infinite time interval can be obtained [CLIO]. 

Open Problem: Are polynomial bounds imposed on g(s) and go{s) and strong 
coercivity condition imposed on go necessary for obtaining weak solutions? 
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3.8 KirhhofF-Boussinesq equation with boundary 
damping - model (2.21) 

Rotational case a > 

In the case when rotational inertia are retained in the model, a > 0, the nonlin- 
ear terms are subcritical and the well-posedness theory can be carried out along 
the same lines as in the case of von Karman equations. For instance, when 
considered (2.21) with g — 0, G ~ and the following boundary damping 

w = 0, Aw = —goi-g-Wt) on S, (3.36) 

where gois) is continuous and monotone, full Hadamard local well-posedness 
can be established. 

Theorem 3.22 Under the assumption 3.19 with g ~ and a > for all initial 
data 

w{o) ^woe H^{n)r\H^in),wt{0) = wi e H^in) 

there exist a unique, local ( in time ) generalized ( semigroup ) solution w of finite 
energy, i.e., 

3T>0: weCi[0,T];H\n)nH^{n)), wt e C{[0,T]; H^in)). 
Moreover 

• //, in addition, the damping go{s) are of some polynomial growth and 
{goix) — goiy)){x — y) > c\x — yY for some r > \, then a generalized 
solution becomes also weak solution. Weak solution is also continuously 
dependent on the initial data. 

• The solutions are global under the dissipativity hypothesis in (3.24) 

The proof of this theorem is along the same lines as in the case of von Karman 
equations with boundary damping. 

Non-rotational case: a = 

As we have seen already before, this case is much more delicate even in the 
case of interior damping. The nonlinear source term is supercritical. While 
this difficulty was overcome in the case of linear interior damping, the presence 
of boundary damping brings new set of issues. This is both, at the level of 
uniqueness and continuous dependence. More specifically, in the case 5 = 0: 

• Existence of finite energy solutions 

w e L.^i[0,T];H'\n) n H^m, wt G L^{[0,T]; L^in)) 
defined variationally and locally in time with initial data 

w{0) = wae H^{n) n H^{n),wt{0) = e L2{n). 
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can be proved by Galerkin method with a boundary damping go{s) of 
some polynomial growth and such that (gaix) — gQ{y)){x — y) > c\x — 
for some ?' > 1. 

• Uniqueness of solutions can be established only in the case of linear damp- 
ing, i.e., for 50 (s) = as. This can be done by adapting Sedenko's method 
[Sed91] as in [Las98, KL02] where full Von Karman system was considered. 
The latter displays similar difficulties when dealing with well-posedness. 

• Continuous dependence on the data can be proved only in a weak topology. 
The difficulty lies in the fact that time reversibility of the flow is lost with 
the presence of boundary damping. This latter property was critical in 
proving energy identity for weak solutions in the case of internal linear 
damping. 

Generalizations- Extensions 

1. Existence and uniqueness of regular solutions in H'^iyi) x H'^{Vl) (when 
a = 0). This can be accomplished by taking advantage of higher topologies 
for the state, hence avoiding the problem of supercriticality. This method 
has been pursued in [CLll] in the case of interior damping. However, the 
arguments are applicable to the case of boundary damping as well. 

2. Other types of boundary damping -such as ocurring in "free" [CLIO] 
boundary conditions can also be considered. In the case rotational case, 
well-posedness results are complete - in line with Theorem 3.22. In the 
non-rotational case, comments presented above apply to this case as well 
(with the same limitations regarding linearity of the damping). There is 
however an additional difficulty resulting from the fact that free boundary 
conditions provide intrinsically nonlinear contribution on the boundary. 
This, however, is of lower order and can be handled as in [CL06b, CLll] 

3. Combination of internal damping and boundary damping can be treated 
by combining the results available for each case separately. 

4. Sources that are nonconservative also can be considered, but these are 
outside the scope of these lectures, see [CLIO]. 

Open Problems: 

1. Hadamard well-posedness with linear boundary damping in the case a = 0. 
Due to the loss of time reversibility, the arguments used for the internal 
damping arc no longer applicable. 

2. Ultimately: uniqueness and Hadamard well-posedness in the presence of 
nonlinear damping go in the case a = and boundary conditions (3.36), 
where go{s) is monotone and -say- of linear growth at infinity. The argu- 
ments used so far for the uniqueness of weak solutions rely critically on 
linearity of the damping. Thus, the case described above appears to be 
completely open. 
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4 General tools for studying attractors 



In this section we describe several approaches to the study of long-time behavior 
of hyperbolic-like systems described above. For the general discussion of long- 
time behavior of systems with dissipation we refer to the monographs [BV92, 
Chu99, Hal88, Lad91, SY02, Tem88]. 

4.1 Basic notions 

By definition a dynamical system is a pair of objects {X, St) consisting of a 
complete metric space X and a family of continuous mappings {St ; t € R-i-} 
of X into itself with the semigroup properties: 

So = I, St+T = St O Sr- 

Wc also assume that y{t) ~ Styo is continuous with respect to t for any i/q €z X . 
Therewith X is called a phase space (or state space) and St is called an evolution 
semigroup (or evolution operator). 

Definition 4.1 Let {X,St) be a dynamical system. 

• A closed set B C X is said to be absorbing for {X, St) iff for any bounded 
set D G X there exists to{D) such that StD C B for all t > tQ{D). 

• (X,St) is said to be (bounded, or ultimately) dissipative iff it possesses 
a bounded absorbing set B. If X is a Banach space, then a value R > is 
said to be a radius of dissipativity of {X, St) if B d {x Cz X : \\x\\x < R}- 

• {X,St) is said to be asymptotically smooth iff for any bounded set D 
such that StD C D for t > Q, there exists a compact set K in the closure 
D of D, such that 

lim dx{StD | if} = 0, (4.1) 
where dx{A\B} = sup^g^ distx(2;, S). 

A set D G X IS said to be forward (or positively) invariant iff StD C D for 
all t > Q. It is backward (or negatively) invariant iff StD D D for all i > 0. 
The set D is said to be invariant iff it is both forward and backward invariant; 
that is, StD = D for all t > 0. 

Let D cX. The set 

7^, ^ U SrD 

T>t 

is called the tail (from the moment t) of the trajectories emanating from D. It 
is clear that 7^ = 7stD- 

If Z3 = {v} is a single point set, then 7+ := 7^ is said to be a positive 
semitrajectory (or semiorbit) emanating from v. A continuous curve 7 = {u{t) : 
i G M} in X is said to be a full trajectory iff Stu{T) = u{t + r) for any t e M 
and t > 0. Because St is not necessarily an invertible operator, a full trajectory 
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may not exist. Semitrajectories are forward invariant sets. Full trajectories are 
invariant sets. 

To describe the asymptotic behavior we use the concept of an w-limit set. 
The set 

t>0 t>QT>t 

is called the uj-limit set of the trajectories emanating from D (the bar over a 
set means the closure). It is equivalent to saying that x e ^{D) if and only if 
there exist sequences t„ +00 and Xn €z D such that St^Xn — > a; as n — > 00. 
It is clear that w-limit sets (if they exist) are forward invariant. 

4.2 Criteria for Asymptotic Smoothness 

The following assertion is a generalization Ceron-Lopes criteria (see [HalSS] and 
the references therein). 

Theorem 4.2 Let (X,St) be a dynamical system on a Banach space X. As- 
sume that for any bounded positively invariant set B in X there exist T > 0, 
a continuous nondecreasing function g : M+ H- R+ , and a pseudometric g'^ on 
C{0,T;X) such that 

(i) 5(0) ^ 0; g{s) <s, s>Q. 

(a) The pseudometric gj^ is precompact (with respect to the norm of X) in the 
following sense. Any sequence {x„} C B has a subsequence {a;„j,} such 
that the sequence {yk} C C(0,T; X) of elements ^^(t) — SrXn^. is Cauchy 
with respect to g]^. 

(Hi) The following estimate 

WSryi - STy2\\ < 9 (||yi - yi\\ + &{{Sry\}ASryiY)) 

holds for every j/i,y2 G B, where we denote by {Sryi} the element in the 
space C(0,T;X) given by function yi{T) = Sryi- 

Then, {X,St) is asymptotically smooth dynamical system. 

Note that a precompact pseudometric is evaluated on trajectories St-, rather 
than on initial conditions (as in the classical treatments, see, e.g., [HalSS]). 
This fact becomes quite useful when applying the criterion to hyperbolic-like 
dynamics. 

Proof. We refer to [CLOSa] . The main ingredient of the proof is the relation 

aiSrB) < g{a{B)), (4.2) 

where ol{B) is the Kuratowski a-measure of noncompactness which is defined 
by the formula 

a{B) = mi{5 : B has a finite cover of diameter < S} (4.3) 
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for every bounded set B of X. For properties of this metric characteristic we 
refer to [Hal88] or [SY02, Lemma 22.2]. 

The property in (4.2) imphes that for every bounded forward invariant set 
B there exists T > such that aiSrB) < a{B) provided a{B) > 0. If for 
some fixed T > this property holds for every bounded set B such that StB 
is also bounded, then, by the definition (see [Hal88]), St is a conditional a- 
condensing mapping. It is known [Hal88] that conditional a-condensing map- 
pings are asymptotically smooth. Therefore Theorem 4.2 can be considered as 
a generalization of the results presented in [Hal88] . ■ 

The above criterion is rather general, however it requires "compactness" of 
the sources in the equation. There are two other criteria that avoid such a 
requirement of a priori compactness. These are: 

• compensated compactness criterion -an idea introduced in [Kha06] and 
later expanded in [CL08a]; 

• J. Ball's energy method, see [Bal04] and [MRW98]. 
The corresponding results are presented below. 



Compensated compactness method 

Theorem 4.3 Let {X, St) be a dynamical system on a complete metric space 
X endowed with a metric d. Assume that for any bounded positively invariant 
set B in X and for any e > there exists T = T(e, B) such that 

d{STyi, STy2) < e + ^'€,B,T(yi, 2/2), yi e B, (4.4) 

where t.B,T{yi,y2) is a functional defined on B x B such that 

liminf liminf T(yr, , y^ri) =0 (4.5) 

for every sequence {y-n} from B. Then {X,St) is an asymptotically smooth 
dynamical system. 

Note that in a "compact" situation, i.e., when the functional ^ is sequentially 
compact, the condition (4.5) is automatically satisfied. The above criterion 
applies in the case of critical nonlinearities. 

Proof. The properties in (4.4) and (4.5) makes it possible to prove that 
a{StB) — >■ as i — )■ -l-oo, where a{B) is the Kuratowski a-mcasurc defined 
by (4.3). The latter property implies the conclusion desired. For details we 
refer [CL08a] or [CLIO]. ■ 
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John Ball's "energy" method 

This second method appHes even in the case of supercritical nonlinear terms, 
but there are other requirements which restrict applicability of this method. 
The main idea behind Ball's method is to construct an appropriate energy type 
functional which can be then decomposed into exponentially decaying part and 
compact part. While the idea of decomposition of semigroup into uniformly 
stable and compact part is behind almost all criteria leading to asymptotic 
smoothness, the "energy" method described below postulates such a decom- 
position on functionals rather than operators (semigroups). This fact has far 
reaching consequences and allows application of the method in supercritical 
situations. 

We follow presentation of the method given in [MRW98] . Another exposition 
of this method in the case of the damped wave equation can be found in [Bal04] . 

Theorem 4.4 Let St be an evolution continuous semigroup of weakly and stron- 
gly continuous operators. Assume that there exist functionals $, ^! , L, K on 
phase space such that the following equality 



holds for any u £ X. The functionals <f>, L, K are assumed to have the 
following properties: 

• <t> : X ^ K.^ is continuous, hounded and if {Uj}j is hounded in X, 
tj +00, StjUj U weakly in X, and limswp^_^^ ^{StjUj) < ^{U), 
then St-Uj — > U strongly in X. 

• ^' : X I— > M is 'asymptotically weakly continuous' in the sense that if{Uj}j 
is hounded in X, tj — > +00, St^Uj U weakly in X, then "^{StjUj) — > 



m K : X I— > R is 'asymptotically weakly continuous' in the sense that if 
{Uj}j is hounded in X, tj — >■ +00, St^ Uj U weakly in X , then K{SsU) G 
ii(0,t) and 



• L is 'asymptotically weakly lower semicontinuous' in the sense that if 




(4.6) 




{Uj}j is hounded in X , tj 
L{SsU) G ii(0,t) and 



+00, StjUj U weakly in X, then 
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Then St is asymptotically smooth. 

Remark 4.1 1. The decomposition of the functionals in (4.6) depends on 
vahdity of energy identity satisfied for weak solutions. This, alone, is a 
severe condition which may be difficult to verify (in contrast to energy 
inequality). In addition, the proof of equality in (4.6) hides behind the 
fact that the damping in the equation is very structured-typically linear. 

2. The functional (typically convex) $ plays role of the energy of linearized 
system - a good topological measure for the solution. In uniformly convex 
spaces X the assumptions postulated by ^ are automatically satisfied. In- 
deed, this results from the fact that weak convergence and the convergence 
of the norms to the same element imply strong convergence. 

3. The hypotheses required from K and L represent some compactness 
property of part of the nonlinear energy describing the system. This is 
often the case even for supercritical nonlinearities. In fact, these terms 
allow to deal with non-compact sources in the equation provided that 
the corresponding nonlinear part of the energy is sequentially compact. 
Typical application involves 3D wave equation, where sources |/(s)| < 
C(l -I- \s\P) with p < 5 can be handled (in view of compactness of the 
embedding H^{n) C Lg{n), q <6). 

4.3 Global attractors 

The main objects arising in the analysis of long-time behavior of infinite-dimen- 
sional dissipative dynamical systems are attractors. Their study allows us to 
answer a number of fundamental questions on the properties of limit regimes 
that can arise in the systems under consideration. At present, there are sev- 
eral general approaches and methods that allow us to prove the existence and 
finite-dimensionality of global attractors for a large class of dynamical systems 
generated by nonlinear partial differential equations (see, e.g., [BV92, Chu99, 
Hal88, LadQl, Tem88] and the references listed therein). 

Definition 4.5 A bounded closed set A G X is said to be a global attractor of 
the dynamical system {X, St) iff the following properties hold. 

(i) A is an invariant set; that is, St A = A for t > 0. 

(ii) A is uniformly attracting; that is, for all bounded set D <Z X 

lim dx{StD\A} ^0, 

where dx{A\B} = sup^^^ distx(a;, _B) is the Hausdorff semidistance. 

It turns out that dissipativity property along with asymptotic smoothness 
imply an existence of global attractor. The corresponding result is standard by 
now and reported below (see [IIal88] and also [BV92, LadQl, TemSS]). 



37 



Theorem 4.6 Any dissipative asymptotically smooth dynamical system {X, St) 
in a Banach space X possesses a unique compact global attractor A. This attrac- 
tor is a connected set and can be described as a set of all hounded full trajectories. 
Moreover A — (^{B) for any hounded absorbing set B of {X, St). 

In the case when the dynamical system has special property-referred to as 
gradient system- dissipativity property is not needed (in the explicit form) in 
order to prove existence of a global attractor. This is a very handy property 
particularly when the proof of dissipativity is technically involved. 

Gradient systems 

The study of the structure of the global attractors is an important problem from 
the point of view of applications. There are no universal approaches solving this 
problem. It is well known that even in finite-dimensional cases an attractor can 
possess extremely complicated structure. However, some sets that belong to 
the attractor can be easily pointed out. For example, every stationary point 
{Stx = X for all t > Q) belongs to the attractor of the system. One can shows 
that any bounded full trajectory also lies in the global attractor. 
Wc begin with the following definition. 

Definition 4.7 Let J\f be the set of stationary points of the dynamical system 
iX,St): 

JV = {veX : Stv^v for all t > 0} . 

We define the unstable manifold Ai'^{Af) emanating from the set Af as a set 
of all y € X such that there exists a full trajectory j — {u{t) : i e M} with the 
properties 

u{Q) ^ y and Um distx{u{t),jV) = 0. 

Now we introduce the notions of Lyapunov functions and gradient systems (see, 
e.g., [BV92, Chu99, Hal88, Lad91, TcmSS] and the references therein). 

Definition 4.8 Let Y Q X be a forward invariant set of a dynamical system 
{X,St). 

• The continuous functional $(y) defined on Y is said to be the Lyapunov 
function for the dynamical system {X,St) on Y iff the function t i— >■ 
^(Sty) is a nonincreasing function for any y Cz Y . 

• The Lyapunov function $(y) is said to be strict on Y iff the equation 
$(5(2/) = ^(y) for all t > and for some y € Y implies that Sty ~ y for 
all t > 0; that is, y is a stationary point of {X, St). 

• The dynamical system {X, St) is said to be gradient iff there exists a 
strict Lyapunov function for {X, St) on the whole phase space X . 
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A connection between gradient systems and existence of compact attractors is 
given below. The main result stating existence and properties of attractors for 
gradient systems is the following theorem (for the proof we refer to [CL08a, 
CLIO] and the references therein; see also [Rau02, Theorem 4.6] for a similar 
assertion) . 

Theorem 4.9 Assume that {X, St) is a gradient asymptotically smooth dynam- 
ical system. Assume its Lyapunov function $(a;) is bounded from above on any 
bounded subset of X and the set ^r~{x : 3>(a;) < R} is bounded for every R. 
If the set J\f of stationary points of {X, St) is bounded, then {X, St) possesses a 
compact global attractor A = M'^{M). Moreover, 

• The global attractor A consists of full trajectories 7 = {u{t) : t G M} such 
that 

lim Aisi X {u{t), N) and lim distx(u(t),7V) = 0. 

• For any x Cz X we have 

lim dist X (Stx, Af) ^ 0, 

that is, any trajectory stabilizes to the set M of stationary points. In 
particular, this means that the global minimal attractor A^un coincides 
with the set of the stationary points, A^m^ = J\f. 

• If N = {zi, . . . , Zn} is a finite set, then A = ^'^^iM'^ {zi) , where M'^{zi) 
is the unstable manifold of the stationary point z^, and also 

(i) The global attractor A consists of full trajectories 7 = {u{t) : t € M} 
connecting pairs of stationary points: any u Cz A belongs to some full 
trajectory 7 and for any 7 C A there exists a pair {z, z*} C Af such 
that 

u{t) — >■ z as t — > —00 and u{t) — > z* as t ^ +00. 
(a) For any v X there exists a stationary point z such that Stv z as 

t — > +CX). 

Dimension of global attractor 

Finite-dimensionality is an important property of global attractors that can be 
established for many dynamical systems, including those arising in significant 
applications. There are several approaches that provide effective estimates for 
the dimension of attractors of dynamical systems generated by PDEs (see, e.g., 
[BV92, LadQl, Tem88]). Here we present an approach that does not require 
C^-smoothness of the evolutionary operator (as in [BV92, TemSS]). The reason 
for this focus is that dynamical systems of hyperbolic-like nature do not display 
smoothing effects, unlike parabolic equations. Therefore, the smoothness 
of the flows is most often beyond question, particularly in problems with a 
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nonlinear dissipation. Instead, we present a method which can be apphed to 
more general locally Lipschitz flows. This method generalizes Ladyzhenskaya's 
theorem (see, e.g., [Lad91]) on finite dimension of invariant sets. We also refer to 
[Pra02] for a closely related approach based on some kind of squeezing property. 
However, we wish to point out that the estimates of the dimension based on the 
theorem below usually tend to be conservative. 

Definition 4.10 Let M be a compact set in a metric space X . 

• The fractal (box-counting) dimension dim^Af of M is defined by 



where n{M, e) is the minimal number of closed balls of the radius e which 
cover the set M . 

Wc can also consider the Hausdorff dimension dim/f to describe complexity and 
embeddings properties of compact sets. We do not give a formal definition of 
this dimension characteristic (see, e.g., [Fal90] for some details and references) 
and we only note that (i) the Hausdorff dimension does not exceed (but is not 
equal, in general) the fractal one; (ii) fractal dimension is more convenient in 
calculations. 

The following result which generalizes [Lad91] was proved in [CL04a], see 
also [CL08a]. 

Theorem 4.11 Let H be a separable Banach space and M be a bounded closed 
set in H . Assume that there exists a mapping V : AI i-^ H such that 

(i) M C VM. 

(ii) V is Lipschitz on M; that is, there exists L > such that 



(Hi) There exist compact seminorms ni{x) and n2{x) on H such that 

\\Vvi - Vvi\\ < ri\\vi - V2\\ + K ■ [ni{vi - V2) + n2{Vvi - VV2)] 

for any Vi,V2 € M , where < 77 < 1 and K > are constants (a seminorm 
n{x) on H is said to be compact ifjn{xm) — ^ for any sequence {xm} C H 
such that Xm weakly in H). 

Then M is a compact set in H of a finite fractal dimension. Moreover, if H 
is a Hilbert space and the seminorms ui and n2 have the form ni{v) ~ \\Piv\\, 
i — 1,2, where Pi and P2 are finite- dimensional orthoprojectors, then 



dmif M = limsup 



lnn(M,e) 
ln(l/£) 



\\Vvi~Vv2\\ < L\\vi-V2\\, vi,V2 e M. 
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Remark 4.2 We note that under the hypotheses of Theorem 4-11 the mapping 
V possesses the property (see Lemma 2.18 in [CL08a]) 

a{VB) < r]a{B) for any B C M, 

where a{B) is the Kuratowski a-measure given by (4.3). This means that V is 
a-contraction on M (in the terminology of [Hal88]). The note that the latter 
property is not sufficient for finite- dimensionality of the set M . Indeed, let 

X = l2 = ^ = {xi]X2] ■■■) : ^2:^ < oo| 

and 

M ^ [x = {xi;x2; ■ ■ ■) &h ■ \xi\<i^^, i = l,2,...} 
We define a mapping V in X by the formula 

[Vx]^ = f,ix^,), 1 = 1,2,... 

where fi{s) = s for \s\ < i^"^ , fi{s) = i^^ for s > i'"^ , and fi{s) = — for 
s < —i~^. One can see that V is globally Lipschitz in X and VX = M = VM. 
Since M is a compact set, the mapping is a-contraction (with rj = 0). On the 
other hand it is clear that dim / M = 00 . 

We also not that this example means that the statement of Theorem 2.8.1 in 
[Hal88] is not true and thus it requires some additional hypotheses concerning 
the mapping. 

In what follows we shall present a unified criterion which allows to prove 
both finite-dimensionality and smoothness of attractors. This criterion involves 
a special class of systems that we call quasi-stable. 

4.4 Quasi-stable systems 

In this section following the presentation given in [CLIO, Section 7.9] we intro- 
duce a class of dissipative dynamical systems which display rather special long 
time behavior dynamic properties. This class will be referred to as quasi-stable 
systems which are defined via quasi- stability inequality given in Definition 4.13. 
It turns out that this is quite large class of systems naturally occurring in non- 
linear PDE's of hyperbolic type of second order in time possibly interacting with 
parabolic equation. The interest in this class of systems stems from the fact that 
quasi- stability inequality almost automatically implies -in one shot- number of 
desirable properties such as asymptotic smoothness, finite dimensionality 
of attractors, regularity of attractors, exponential attraction etc. In 
what follows we shall provide brief introduction to the theory of quasi-stable 
systems. 

We begin with the following assumption. 
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Assumption 4.12 Let X, Y , and Z he reflexive Banach spaces; X is compactly 
embedded in Y . We endow the space H = X x Y x Z with the norm 

IvIh = \uo\x + \ui\y + l^lli y = {uq;ui;0q). 

The trivial case Z ~ {0} is allowed. We assume that {H,St) is a dynamical 
system onH = XxYxZ with the evolution operator of the form 

Sty = {u{t);ut{t);e{t)), y ^ {uo;ui;0o) £ H, (4.7) 

where the functions u(t) and 9(t) possess the properties 

ueC{R+,x)r\C\R+,Y), eeC{R+,z). 

The structure of the phase space H and the evolution operator St in Assump- 
tion 4.12 is motivated by the study of the system generated by equation of 
the second order in time in A x 1^ possibly interacting with some first-order 
evolution equation in space Z . This type of interaction arises in modeling of 
thermoelastic plates and structural acoustic sytems. 

Definition 4.13 A dynamical system of the form (4.7) is said to be stable mod- 
ulo compact terms (quasi-stable, for short) on a set B d H if there exist a 
compact seminorm fJ.xi') on the space X and nonnegative scalar functions a(t), 
b{t), and c{t) on M-|_ such that (i) a{t) and c{t) are locally bounded on [0,oo), 
(ii) bit) G Li(]R_|-) possesses the property limt^oo 6(i) = 0, and (Hi) for every 
2/1 J 2/2 G /S and t > Q the following relations 

\Styi-Sty2\l<a{t)-\yi-y2\l (4.8) 

and 

\Styi-Sty2?H<b{t)-\yi-y2^H + <t)- sup [^,x{u\s) - u'is))]' (4.9) 

0<s<t 

hold. Here we denote Styi = {u^{t);ul{t);0'^{t)), i = 1,2. 

Remark 4.3 The definition of quasi- stability is rather natural from the point 
of view of long-time behavior. It pertains to decomposition of the flow into 
exponentially stable and compact part. This represents some sort of analogy 
with the "splitting" method [BV92, Tem88], however, the decomposition refers 
to the difference of two trajectories, rather than a single trajectory. In addition, 
the quadratic dependence with respect to the semi-norm in (4.9) is essential in 
the definition. The relation (4.9) is called a stabilizability estimate (or quasi- 
stability inequality) and, in the context of long-time dynamics, was originally 
introduced in [CL04a] (see also [CEL04, CL04b] and the discussion in [CL08a] 
and [CLIO]). To obtain such an estimate proves fairly technical (in critical 
problems) and requires rather subtle PDE tools to prove it. Illustrations of the 
method are given in [CLIO] for some abstract models and also for a variety of von 
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Karman models. We also refer to [BC08, BCL07, CEL04, CL04a, CL06a, CL07a, 
CL08a, CLT08, CLT09, Nab09] for similar considerations for other models. 

The notion of quasi-stability introduced in Definition 4.13 requires a special 
structure of the semiflow and a special type of a (quasi-stability) inequality. 
However, the idea behind this notion can be applied in many other cases (see, 
e.g., [CL04b, CLOSa, CLll] and also [CLIO]. Systems with delay/memory terms 
can be also included in this framework (see, e.g., [Fas07, Fas09, Pot09, Ryz05] 
and also [CLIO]). The same idea was recently applied in [CKIO] for analysis of 
long-time dynamics in a parabolic- type model (see below in Section 6). 

Remark 4.4 In order to write down a more explicit form of quasi-stability 
inequality let us consider dynamical system (say of gradient form) generated 
by some second order evolution equation in the space H = X x Y with an 
associated energy functional -say Ey{t) and the damping integral denoted by 
Dy{t), so that the energy identity reads 

Ey{t) + Dy{T)dT = Ey{s), S<t. 

In this case a sufficient condition for quasi-stability is the following [CEL04]sia- 
bilizability inequality: there exists a parameter T > 0, such that the difference 
of two any trajectories z{t) = yi{t) — y2{t) satisfies the relation 

EoAT)+ [ EoAr)dT<CT [ D,{t) + CtLOT,, (4.10) 
Jo Jo 

where Eg^zit) denotes a positive part of the energy, Dz is the damping func- 
tional for the difference z and LOT^ is a lower order term of the form LOT^ = 
sup^g(-Q ||z(t)|||^^ , where H C Hi with a compact embedding. If one is 
already equipped with an existence of compact attractor, then stabilizabil- 
ity inequality that is required is of milder form: there exist two parameters 
—oo<s<T such that 

EoAT) + J^ EoAr)dT < Cs,T Dz{T)dT + C.^tLOT,. (4.11) 

Here LOT^ denotes lower order terms measured by 

LOT, = sup ||z(t)||^^, 

re(s,T) 

and parameters s,r may depend on the specific trajectories 2/1,2/2 • 

The advantage of the above formulation is that these are inequalities that are 
obtainable by multipliers when studying stabilization and controllability (see, 
e.g., the monographs [Las02, LTOO] and the references therein). The stabiliz- 
ability inequality in (4.10) provides a direct link between controllability and 
long time behavior. Indeed, when considering stabilization to zero equilibria. 
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the lower order terms LOTz are dismissed by applying compactness-uniqueness 
argument. Then the inequality applied to a single trajectory yi = y 

Ey{T) < Co,T r Dy{T) = Co.T[Ey{0) - Ey{T)] (4.12) 
Jo 

implies, by standard evolution method, exponential decay to zero equilibrium. 

Similarly, in the case of controllability the observability inequality required 
[Las02] is precisely (4.12) with Dy{t) denoting the quantities observed (for in- 
stance, velocity in the spatial domain or velocity on the subdomain). In cither 
case, the crux of the matter is to establish such inequality. 

In what follows our aim is to show that quasi-stable systems have far reach- 
ing consequences and enjoy many nice properties that include (i) existence 
of global attractors that is both finite-dimensional and smooth, (ii) 
exponential attractors, and so on. 

Asymptotic smoothness 

Proposition 4.14 Let Assumption ^.12 be in force. Assume that the dynamical 
system {H,St) is quasi-stable on every bounded forward invariant set B in H. 
Then, [H, St) is asymptotically smooth. 

Proof. Let 

X = Closure {v<eX : \v\^ = ^ix{v) + \v\y < oo] . (4.13) 

One can sec that X is compactly embedded in the Banach space X. We have 
that the space 

M^i,2(0,T;X,y) = {/ e Loo(0,T;X) : /' € MQ^T-Y)} 

is compactly embedded in C{Q,T;X). This implies that the pseudometric q*^ 
in C(0, t; H) defined by the formula 

QB{{Sryi}, {S'r?/2}) = c{t) SUp ^lx{u^{t) - {t)) 

re[o,t] 

is precompact (with respect to H). Here we denote by {Sryi} the element from 
C{{),t-H) given by function y,{T) = Sry^ = {u'{t);ul{t);e'{t)). By (4.9) g'g 
satisfies the hypotheses of Theorem 4.2. This implies the result. ■ 

Corollary 4.15 If the system [H, St) is dissipative and satisfies the hypotheses 
of Proposition 4.14, then it possesses a compact global attractor. 

Proof. By Proposition 4.14 the system {H, St) is asymptotically smooth. Thus 
the result follows from Theorem 4.6. I 
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Finite dimension of global attractors 

We start with the following general assertion. 

Theorem 4.16 Let Assumption 4.12 be valid. Assume that the dynamical sys- 
tem {H, St) possesses a compact global attractor A and is quasi-stable on A (see 
Definition 4.13J. Then the attractor A has a finite fractal dimension dimj A 
(this also implies the finiteness of the Hausdorff dimension) . 

Proof. The idea of the proof is based on the method of "short" trajectories 
(see, e.g., [MN96, MP02] and the references therein and also [CL08a]). 

We apply Theorem 4.11 in the space Ht = HxWi{Q, T) with an appropriate 
T. Here 

W,{0,T) = |z e L2(0,T;X) : |z|^,(o,t) =^ {W)\\ + NtWIy) < ooj . 

(4.14) 

The norm in Ht is given by 

llC^ilk = + l^lk(o,T)> U ={y-z), y={uQ;ui-eQ). (4.15) 
Let Hi = (uq; u\] 9q), i = 1,2, be two elements from the attractor A. We denote 

Sty, = iu\t);ul{t);9\t)), t>0,i^l,2, 
and Z{t) = Styi - Sty2 = {z{t);zt{t);£,{t)), where 

{z{t); zt{t);m) ^ {u\t) - u\t); u\{t) - ul{t);e\t) - e\t)). 
Integrating (4.9) from T to 2T with respect to t, we obtain that 

/ \Styi- Sty2\%dt<bT\yi-y2\]i + CT sup [^ix{z{s))f , (4.16) 

JT 0<s<2T 

where 

_ p2T p2T 

bT= I b{t)dt and ct ^ c{t)dt. 

JT JT 

It also follows from (4.9) that 

\STyi - STy2\]i < b{T) ■ \yi - y2\'\j + c{T) ■ sup [^ix{z{s))f 

0<s<T 

and combining with (4.16) yields 



\STyi-STy2?H+ I \Styi - Sty2\Hdt < bT\yi - y2\H + CT sup [fix{z{s))]'^ , 

JT 0<s<2T 

(4.17) 



where 

r.2T r2T 



bT = b{T) + / b{t)dt and ct = c{T) + / c{t)dt. (4.18) 

JT JT 
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Let A be the global attractor. Consider in the space Ht the set 

At := {U = {u{0); ut{0); 0(0); t G [0, T]) : (u(0); ut{0); 6(0)) e A} , 

where u{t) is the first component of Sty{0) with ?;(0) = (m(0); Ut(0); 0(0)), and 
define operator V : At i-> Ht by the formula 

V : {u{0);ut{0);9{0);u{t)) ^ {STy{0);u{T + t)). 

It is clear that V is Lipschitz on At and V^Ay = At- 

Because the space X given by (4.13) possesses the properties 

X C X C Y and X C X is compact, 

contradiction argument yields 

[^ix{u)f <e\u\l + CMY for any e > 0. (4.19) 

Therefore it follows from (4.8) that 

sup [fixiz{s))f <—\yi-y2\H + C{aT,bT,CT) sup |z(s)|y, 

0<s<2T Ct 0<s<2T 

where ut = supQ<g<2T Q^(s) and hT,CT given by (4.18). Consequently, from 
(4.17) we obtain 

||y;7i - VU2\\ht < Vt\\Ui - C/zHh^ + Kt ■ [nriUi - U2) + nT{VUi ~ VU2)). 

for any J7i, C/2 G ^t, where ify > is a constant (depending on aT,bT,CT, the 
embedding properties of X into Y , the seminorm /ix), and 

ri^ = bT = h{T) + / b{t)dt. (4.20) 

The seminorm tit has the form nT{U) := supq<<,<7^ |u(s)|y. Because Wi(0,T) 
is compactly embedded into C{Q,T\Y)^ nT{U) is a compact seminorm on Ht 
and we can choose T > 1 such that tjt < 1- We also have from (4.8) that 

\\VUi-VU2\\h^ <LT\\Ui-U2\\nT for Ui,U2 e At, 

where = a(T) + J^^ a{t)dt. Therefore we can apply Theorem 4.11 which 
implies that At is a compact set in Ht of finite fractal dimension. 
Let V : Ht H he the operator defined by the formula 

A — VAt and V is Lipshitz continuous, thus dim^A < dimj'^ At < 00. Here 
dimJT stands for the fractal dimension of a set in the space W. This concludes 
the proof of Theorem 4.16. ■ 
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Remark 4.5 By [CL08a] the dimension diniy A of the attractor admits the 
estimate 

dimf A < 

Here mo{R) is the maximal number of pairs {xi; jji) in Ht x Ht possessing the 
properties 

W^AIht + WV'-W'ht - ' nT{xt-Xj) + nT{yi-yj)>l, i ^ j- 

It is clear that mQ{R) can be estimated by the maximal number of pairs (x^; y,) 
in Wi{0,T) X Wi{0,T) possessing the properties j,, + < R"^ 

and nxixi — Xj) + nT{yi — yj) > 1 for all i ^ j. Thus the bound in (4.21) 
depends on the functions a, 6, c and scminorm fix in Definition 4.13 and also 
on the embedding properties of X into Y. 

Regularity of trajectories from the attractor 

In this section we show how stabilizability estimates can be used in order to 
obtain additional regularity of trajectories lying on the global attractor. The 
theorem below provides regularity for time derivatives. The needed "space" 
regularity follows from the analysis of the respective PDE. It typically involves 
application of elliptic theory (sec the corresponding results in [CLIO, Chap- 
ter 9]). 

Theorem 4.17 Let Assumption 4.12 be valid. Assume that the dynamical 
system {H,St) possesses a compact global attractor A and is quasi-stable on 
the attractor A. Moreover, we assume that (4.9) holds with the function c{t) 
possessing the property Coo = sup^gjj^ c(t) < oo. Then any full trajectory 
{{u{t);ut(t);9(t)) : t £ R} that belongs to the global attractor enjoys the follow- 
ing regularity properties, 

e Loo(IR;X)nC(M;y), uu e LooiR;Y), eteL^{R;Z) (4.22) 

Moreover, there exists R > such that 

\ut{t)\l + \utt{t)\l^ + \9tit)\l < r:', tern, (4.23) 

where R depends on the constant Coa, on the scminorm fix in Definition 4.13, 
and also on the embedding properties of X into Y. 

Proof. It follows from (4.9) that for any two trajectories 

7 = {Uit) = {u{t);ut{t);d{t)) : teM}, 

7* = {u*it)^iu*ity,u:ity,0*it)) -.teR} 

from the global attractor wc have that 

\Z{t)\l < b{t - s)\Z{s)\l + c{t - s) sup ifixizir))]' (4.24) 

s<r<t 



In- 



In niQ 



AKt{1 



(4.21) 
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for all s<t, s,t eR, where Z{t) = U*{t) ~ Uit) and z{t) = u*{t) - u{t). In the 
limit s — oo relation (4.24) gives us that 

\Z{t)\j, < sup [fix{z{T))]^ 

— oo<r<i 

for every i G M and for every couple of trajectories 7 and 7*. Using relation 
(4.19) we can conclude that 

sup \ZiT)\l<C sup |z(r)|^, (4.25) 

— oo<r<t —oo<T<t 

for every t gM. and for every couple of trajectories 7 and 7* from the attractor. 

Now we fix the trajectory 7 and for < < 1 we consider the shifted 
trajectory 7* = 7/1 = {y{t + h) : t e M}. Applying (4.25) for this pair of 
trajectories and using the fact that all terms (4.25) are quadratic with respect 
to Z we obtain that 

sup {\u'^{r)\j, + \u1ir)\l. + \e'^ir)\l}<C sup \u\r)\l^, (4.26) 

— oo<r<t "'^^'J'^* 

where u''{t) = h-^ ■ [u{t + h) ~ u{t)] and e''{t) = h-^ ■ [e{t + h) - e{t)]. On the 
attractor we obviously have that 

\u'\t)\Y<Yj^ \ut{T + t)\YdT<C, teM, 
with uniformity in h. Therefore (4.26) implies that 

\u\t)\\ + \u^,{t)\l + |0f (i)|| < C, t e M. 

Passing with the limit on h then yields relations (4.22) and (4.23). ■ 



Fractal exponential attractors 

For quasi-stable systems we can also establish some result pertaining to (gener- 
alized) fractal exponential attractors. 
We first recall the following definition. 

Definition 4.18 A compact set Aoxp <Z X is said to be inertial (or a fractal 
exponential attractorj of the dynamical system (X,St) iff A is a positively 
invariant set of finite fractal dimension (in X) and for every bounded set D d X 
there exist positive constants tjj, Cu andju such that 

dx{StD I Aexp} = sup distx(5tx, ^cxp) < Cd ■ e^''"^*-*"), t>tD. 
xeD 

If the dimension of A is finite in some extended space we call A generalized 
fractal exponential attractor. 



48 



For more details concerning fractal exponential attractors we refer to [FEN94] 
and also to recent survey [MZ08] . We only note that the standard technical tool 
in the construction of fractal exponential attractors is the so-called squeezing 
property which says [MZ08], roughly speaking, that either the higher modes are 
dominated by the lower ones or that the semifiow is contracted exponentially. 
Instead the approach we use is based on the quasi-stability property which says 
that the semifiow is asymptotically contracted up to a homogeneous compact 
additive term. 

Theorem 4.19 Let Assumption 4.12 be valid. Assume that the dynamical sys- 
tem {H,St) is dissipative and quasi-stable on some bounded absorbing set B. 
We also assume that there exists a space H D H such that t t-^ Sty is Holder 
continuous in H for every y G B; that is, there exist < 7 < 1 and Cb.t > 
such that 

\St^y-St,y\H<CB,T\ti~t2\\ t,,t2e[0,Tl yeB. (4.27) 

Then the dynamical system (H,St) possesses a (generalized) fractal exponential 
attractor whose dimension is finite in the space H. 

In contrast with Theorem 3.5 [MZ08] Theorem 4.19 does not assume Holder 
continuity in the phase space. At the expense of this we can guarantee finiteness 
of the dimension in some extended space only. This is why we use the notion of 
generalized exponential attractors. 

Proof. We apply the same idea as in the proof of Theorem 4.16. 

We can assume that absorbing set B is closed and forward invariant (other- 
wise, instead of B, we consider B' ~ Llt>toStB for large enough, which lies in 
B). 

In the space Ht = H x Wi{0,T) equipped with the norm (4.15), where 
Wi{0,T) is given by (4.14), wc consider the set 

Bt ■■= {U = (u(0); ut{0); 0(0); u{t),t e [0, T]) : (m(0); ut(0); ^?(0)) e B} , 

where u{t) is the first component of Sty{0) with y(0) ~ (u(0); ut(0); 6'(0)). As 
above we define the shift operator V : Bt '-^ Ht by the formula 

V : {u{0);ut{0);9{0);u{t)) ^ {STy{0);u{T + t)). 

It is clear that Bt is a closed bounded set in Ht which is forward invariant with 
respect to V. 

It follows from (4.8) that 

\\VUi - VU2\\%^ < (^aiT) + a{t)dt^ \\Ui - iJsUl^^, f/i, U2 G Bt- 

As in the proof of Theorem 4.16 we can obtain that 
\\VUi - VU2\\h^ < Vt\\Ui - U2\\ht + Kt ■ MUi - U2) + nT{VUi - VU2)) 
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for any Ui,U2 S Bt and for some T > 0, where Kt > is a constant, 
nxiU) := supo<5<7- and 77T < 1 is given by (4.20). Therefore by Theo- 

rem 7.4.2 [CLIO] the mapping V possesses a fractal exponential attractor; that 
is, there exists a compact set At C Bt and a number < g < 1 such that 
dim^^ At < 00, VAt C At, and 

snp{distH^iV''U, At) : U gBt} <Cq'', ^==1,2,..., 

for some constant C > 0. In particular, this relation implies that 

sup{distH{SkTy,A) : ueB} <Cq'', fc=l,2,..., (4.28) 

where A is the projection of At of the first components: 

A^{iuiO);ut{Oy,e{0)) eB : (w(0); ut(0); 0(0); m(<), t G [0, T]) € ^t} ■ 

It is clear that A is a compact forward invariant set with respect to St', that is, 
St A C A. Moreover dim^ A < dim^^ At < 00. One can also see that 

A,^p^U{StA : te[0,T]} 

is a compact forward invariant set with respect to St', that is. StA^xp C Aexp- 
Moreover, it follows from (4.27) that 



dim;^ A^^p < c 



1 + dimf A 



< 00. 



We also have from (4.28) and (4.8) that 

sup{distH(S'ty,Aexp) : u e B} < Ce'^^ t>0, 
for some 7 > 0. Thus A^^p is a (generalized) fractal exponential attractor. ■ 

Remark 4.6 Holder continuity (4.27) is needed in order to derive the finite- 

ncss of the fractal dimension dim^^ ^cxp from the finiteness of dim^ A. We 
do not know whether the same holds true without property (4.27) imposed in 
some vicinity of A^xp- This is because ^exp is a uncountable union of (finite- 
dimensional) sets StA. We also emphasize that fractal dimension depends on 
the topology. Indeed, there is an example of a set with finite fractal dimension 
in one space and infinite fractal dimension in another (smaller) space. 

5 Long time behavior for canonical models de- 
scribed in Section 2 

The goal of this section is to show how to apply abstract methods presented in 
the previous section in order to establish long time behavior properties of dy- 
namical systems generated by PDE's described in Section 2. An interesting fea- 
ture is that though all the models considered are of hyperbolic type-without an 
inherent smoothing mechanism, the long time behavior can be made "smooth" 
and characterized by finite dimensional structures. Our plan is to focus on the 
following topics. 
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1. Control to finite dimensional and smooth attractors of nonlinear wave 
equation 

• with interior fully supported dissipation, 

• with boundary or partially localized dissipation. 

2. Control to finite dimensional and smooth attractors of nonlinear plates 
dynamics 

• von Karman plates with a nonlinear fully supported interior feedback 
control, 

• von Karman and Berger plates with boundary and partially localized 
damping, 

• Kirchhoff-Boussinesq plates with a linear interior feedback control. 

For each model considered we follow the plan: (i) state the assumptions and the 
results, (ii) discuss possible generalizations, (iii) state open problems, (iv) pro- 
vide a sketch of the proofs. 

5.1 Wave dynamics 

In this section we consider the existence problem for finite dimensional and 
smooth attractors for semilinear wave equation with critical (feedback) sources 
with both interior and boundary damping. 

We start with interior damping and source model as described in (2.1). 

Interior damping 

In studying long time behavior we find convenient to collect all the assumptions 
required for the source and damping; 

Assumption 5.1 1. The source f £ C^(M) satisfies the dissipativity condi- 
tion in (3.4) and the growth condition 

\f'{s)\<M\s\^ for \s\>l. 

2. The damping g(s) is increasing function, g(0) — and it satisfies the 
following asymptotic growth conditions: 

< mg\s\^ < g{s)s < Mg\sf for \s\ > 1. (5.1) 

3. The damping coefficient a satisfies a{x) > oq > 0,a; € f2 . 
With reference to the model (2.1) we can take 

/(s) = -s^ + cs^ and g{s) = gis + \s\'^s with pi > and c G M. (5.2) 
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Remark 5.1 l.The power 3 associated with the source / and the power 5 
associated with the damping g are often referred as " double critical" parameters 
for the wave equation with n ~ 3. The reason for this is that the maps u — > f{u) 
from H^{n) into L2{^) and u g(u) from H^{il) into H^^{n) arc bounded 
but not compact. 

2. More general forms of space dependent and localized damping coefficient 
a{x) are considered in [CLT08] 

Under the above conditions the wellposedness results of Theorem 3.3 assert 
existence of a continuous semiflow St which defines dynamical system {'H,St) 
with H = Hq{Q) X L2{fl)- We are ready to undertake the study of long-time 
behavior. Our first main result is given below: 

Theorem 5.2 (Compactness) With reference to the dynamics described by 
(2.1) and zero Dirichlet boundary conditions, we assume Assumption 5.1 for 
the source f and the damping g. Then, there exists a global compact attractor 
A d H for the dynamical system {T-L, St). In addition (H, St) is gradient system. 

To describe the structure of the attractor, we introduce the set of stationary 
points of St denoted by J\f, 

Af^{Ven : StV tor &l\t>0}. (5.3) 

Every stationary point W £ Af has the form W ~ {w; 0), where w ~ w{x) solves 
the problem 

Aw + f{w) = Omn and w = in T. (5.4) 

Under the Assumption 5.1, standard elliptic theory implies that every stationary 
solution satisfies w £ H'^iyi) and the set of all stationary solutions is bounded 
in H^{Vl). In fact, more regularity of stationary solutions can be shown, but 
the above is sufficient for the analysis to follow. 

The next result is a consequence of general properties of gradient systems 
(see, e.g., [BV92, Hal88] and also Section 4.3 above) and asserts that the at- 
tractor A coincides with this unstable manifold. 

Theorem 5.3 (Structure) Under the assumptions of Theorem 5.2 we have 
that 

• A^ i\/"(A/'); 

• limt^+oodist«(S'tVK,7V) ^QforanyW£H = H^{n) x L2{n); 

• if (5.4) has a finite number of solutions^ , then for any V £% there exists 
a stationary point Z = {z,Q) £ M such that StV Z in % as t ^ +00. 

Our second main result read refers to the dimensionality of attractor. 

Theorem 5.4 (Finite Dimensionality) Let f £ C2(R), g £ C^{M) and As- 
sumption 5.1 be in force. In addition, we assume that for all s £ M; 

-"^We note that the property that (5.4) has finitely many solutions is generic. 
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1- \f"{s)\ < C\s\, 

2. 0< m. < g'{s) < C[l + sg{s)]^^^. 

Then the fractal dimension of the global attractor A of the dynamical system 
{T-L, St) is finite. 

With reference to the model (2.1) we still can take / and g as in (5.2), but with 

The proof of Theorem 5.4 follows from quasi-stability property satisfied for 
the model (see Definition 4.13). The exactly same lemma allows to obtain "for 
free" the following regularity result for elements from the attractor. 

Theorem 5.5 (Regularity) Under the assumptions of Theorem 5.4 the at- 
tractor is a bounded set in 

V = {{u,v) e H^(n) X H^{n) \ ^Au + g{v) e L2in), }, 

i.e. there exist constants c; > such that^ 



A C |(u,v) G H\n) X H^n) 



\\u\\i.n + \\v\\i,n < ci; 
\\Au-g{v)\\ < 02, 



A natural question that can be asked in this context is that of existence of 
strong (i.e., corresponding to topology of strong solutions) attractors. Though 
this latter property is technically related to smoothness of attractors, however 
the corresponding result does not follow from Theorem 5.5, unless the damping 
is subcritical. Additional analysis is needed for that. A first question to address 
in this direction is the existence of attractors for strong solutions. This is to 
say that when restricting solutions to regular initial data, the corresponding 
trajectories converge asymptotically to an attractor Ai which typically may 
be smaller than A the latter corresponds to weak or generalized solutions. A 
first step toward this goal is to show dissipativity of strong solutions (which 
of course does not follow at once from dissipativity of weak solutions, unless a 
problem is linear). Fortunately, dsissipativity of strong solution is, again, direct 
consequence of quasi-stability. Thus, we have this property for "free". 

Theorem 5.6 (Dissipativity of strong solutions) Let the assumptions of 
Theorem 5.4 be in force. Then there exists a number Rq > such that for 
any R > we can find t^ > such that 

\\wtt{t)f + \\wt{t)\\l,-i < Rl for all t>tR (5.5) 

for any strong solution w{t) to problem (2.1) with initial data {wq]Wi) from the 
set 

WR = {{wo,wi)eH^{n)xH^{n) : ||wo||i,q + < i?, \\w2\\ < R} , 

(5.6) 

where W2 = g{wi) — Awq. 



See also Theorem 5.7 below which asserts an additional regularity of attractor. 
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Theorem 5.6 refers to strong solutions. The existence of such, is guaranteed 
once initial data are taken from Wr. With additional calculations (using the 
multiplier Awt) one improves the statement and obtains the dissipativity for 
[[^^(Olli n — ^0' t ^ with initial data from H'^^ft). This, in particular, shows 
that strong attractors (for strong solutions) are in H^{^1) nifg(ri) x HQ{n). An 
interesting question is whether the same regularity is enjoyed by weak attractors. 
In other words, whether H'^{n) D i?o(^) x H^i^) ^ -4? The answer to this 
question is given by 

Theorem 5.7 (Regularity) Under the same assumptions as in Theorem 5.4 
the attractor is a bounded set in H^{n) x L2{^1). 

The proof of this theorem given in [KhalO] exploits the so called "backward 
smoothness on attractor" . 

Decay rates to equilibrium: Knowing that every trajectory converges to 
some equilibrium point one would like to know " how fast" ? The answer to this 
question is given by decay rates of asymptotic convergence of solutions to point 
of equilibria. This is the topic we present next. 

We introduce concave, strictly increasing, continuous function : ]R+ i— > ]R_|_ 
which captures the behavior of 17(5) at the origin possessing the properties 

A:o(0) = and s^ + g\s) < koisgis)) for |s| < 1. (5.7) 

Such a function can always be constructed due to the monotonicity of g, see 
[LT93] or Appendix B in [CLIO]. Moreover, on the strength of Assumption 5.1 
there exists a constant c > such that fc(s) = fco(s) + cs is a concave, strictly 
increasing, continuous function k : such that 

fc(0) = and |s|2 < k{sg{s)) for aU s e R. (5.8) 

Given function k we define 

Hois) ^k (^-i^ , Go{s) = c,{I + Ho)-\c2s), (5.9) 
Q{s) ^s-{I + Gor\s), 

where Ci are some positive constants. It is obvious that Q{s) is strictly mono- 
tone. Thus, the differential equation 

^+Q(cr)=0, i>0, CT(0) = CToeM, (5.10) 
at 

admits global, unique solution cr(i) which, moreover, decays asymptotically to 
zero as t ^ 00. With these definitions we are ready to state our next result. 

Theorem 5.8 (Rate of convergence to equilibria) In addition to Assum- 
ption 5.1, we assume that the set of stationary points Af is finite, and every 
equilibrium V — (v;0) is hyperbolic in the sense that the problem 

Aw + f'{v) ■ w = in fl and w — onT (5.11) 
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has no non-trivial solutions. Then, for any Wq = {wq;wi) G Ji, there exists a 
stationary point V = {v;0) such that 

\\S{t)Wo^V\\n<C{a{[tT-'])) , t > 0, (5.12) 

where C and T are positive constants depending on Wq, [a] denotes the integer 
part of a and a(t) satisfies (5.10) with ctq = C{Wo, V) where C{Wo, V) is a con- 
stant depending on ||Wo||-h and \\V\\-h and Q is defined in (5.9). In particular, 
ifg'iO) > 0, then 

\\S{t)Wo-V\\n < Ce-"* 
for some positive constants C and lo depending on Wq. 

Remark 5.2 Since Q{s) is strictly increasing and Q(0) = 0, the rates described 
by the ODE in (5.10) (see, e.g., [LT93] or [CLIO, Appendix B]) decay uniformly 
to zero. The "speed" of decay depends on the behavior of g'{s) at the origin. 
If g'{s) decays to zero polynomially, then by solving the ODE in (5.10), one 
obtains algebraic decay rates for the solutions to cr{t). If, instead (5 = and 
g'{0) > 0, then Q(cr) = aa for some a > and, consequently, the decay rates 
derived from (5.10) are exponential (see [LT93] for details). In the latter case 
using the dissipitavity of strong solutions (see Theorem 5.6) one can prove by 
interpolation that the exponential decay rate holds for strong solutions in a 
stronger (than in H) topology. 

Boundary damping 

In the case when boundary dissipation is the only active mode of damping 
(model in (2.14) and (2.15) with g = 0) one can still prove that the long time 
behavior is essentially the same as in the case of internal damping. The task of 
achieving this goal is much more technical, due to the necessity of propagating 
the damping from the boundary into the interior. This is done by using familiar 
by now "flux" multipliers. However, the resulting analysis and PDE estimates 
are considerably more complicated and often resort to specific unique continua- 
tion property as well as Carleman's estimates (when dealing with critical cases). 
This topic has been considered in [CEL02, CEL04, CLT09]. The corresponding 
results is given below. 

Theorem 5.9 (Boundary Dissipation) With reference to the equations in 
(2.14)and (2.15), where we take g ~ 0, h[w) = ~w, the source f satisfies dissi- 
pativity property imposed in (3.4), where Ai corresponds to the first eigenvalue 
of the Robin problem and such that |/"(s)| < C(l + |s|) for all s G M. The 
boundary damping go is an increasing, differ entiable function, ,9o(0) = and 
satisfies the asymptotic growth condition: 

m < g'ois) < M, \s\ > 1, m, M > 0, (5.13) 

Then 

1. Under the above assumption the statements of Theorem 5.2 and 5.3 hold. 
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2. If, in addition, (7o(0) > then the statements of Theorem 5.4 and also 
Theorem 5.7, Theorem 5.8 (with appropriate modifications) hold true. 

One fundamental difference between boundary damping and interior damping 
is that the restriction of hnear bound at infinity for the boundary damping 
is essential. One dimensional examples [VanOO] disprove validity of uniform 
stability to a single equilibrium when the damping is superlinear at infinity. 

Generalizations 

With reference to interior damping 

1. Wave equation with homogenous Neumann or Robin boundary conditions 
can be considered in the same way. 

2. The same wave model with additional boundary damping and source, 
as described in (2.14) and (2.15). The corresponding analysis is more 
involved. Details are given in [CL07a]. 

3. Possibility of degenerate damping g{s), where the density function a{x) 
describes possibility of degeneracy of the damping. This situation is also 
studied in [CL07a]. 

4. Strong attractors. These are attractors corresponding to strong solutions. 
In the case when boundary source and damping are involved and the 
damping may degenerate, the description of strong attractors is more in- 
volved. These are bounded sets V x iJ^(ri) where H^{n) C V. The details 
are given in [CL07a] . 

5. Under additional assumptions imposed on regularity of the damping g{s) 
and the source /(s) one can prove that the attractors corresponding to 
Theorem 5.7 are infinitely smooth C°°. This can be done by the usual 
boot-strap argument. 

6. Attractors with localized dissipation only. This corresponds to having a{x) 
localized in the layer near boundary. [CLT08] . As in the case of boundary 
damping, requires linearly bounded damping. 

With reference to boundary damping 

1. Attractors with boundary damping which is further localized to a suitable 
portion of the boundary. Requires appropriate geometric conditions, see 
[CLT09] for details. 

2. Relaxed regularity assumptions imposed on the damping function .go(s), 
see [CLT09] again. 
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Open Problems 

1. Higher regularity of attractors with minimal hypotheses on the damping 
and source. 

2. Under which conditions weak attractors coincide with strong attractors in 
the case of boundary or partially localized damping? 

3. Supercritical sources (e.g., /(s) = s^). Ball's method will provide exis- 
tence of attractors with linear damping. However, linear damping is not 
sufficient for uniqueness of solutions. One needs to consider non-unique 
solutions. 

4. Nonmonotone damping. Some particular results in this direction can be 
found in [CL08a]. 

5. Non-gradient structures. For instance models involving first order differ- 
entials as the potential sources. 

Remarks 

1. Theorem 5.2 -Compactness. Difficulty: double critical exponent -rules out 
previous methods. Use Theorem 4.3 after an appropriate rewriting of the 
source which then fits into the structure of the functional This point 
is explained below. 

2. Theorems 5.4 arid Theorem 5.7-finite dimensionality and smoothness. Re- 
lies on proving quasi-stability inequality. The main trick is to provide 
appropriate decomposition of the source that takes advantage of dissipa- 
tivity integral which is Li(R). See [CL07a, CLT09]. 

3. Full statement of Theorem 5.7 - [KhalO]. Smoothness of the trajectories 
on the attractor near — oo. The existence of strong attractors (for strong 
solutions) is de facto consequence of quasi-stability inequality. However, in 
many situation proving additional spatial regularity {H'^{n), for instance) 
is problematic. It is for this task that special considerations related to 
backward smoothness are important. 

4. Decay rates - Theorem 5.8. The difficulty relates to the fact that con- 
vergence to equilibria is a very "unstable" process. It is not possible to 
localize the analysis to the neighborhood of the equilibria. Decay rates 
are derived by using convex analysis and reducing estimates to nonlinear 
ODE'S. See [CL07a]. 

5. Boundary damping - Theorem 5.9. First results are given in [CEL02, 
CEL04] where finite-dimensionality was proved for subcritical sources only. 
Complete proof of Theorem 5.9 (in fact in a more general version) is 
given in [CLT09], where fiux multipliers are used along with Carleman's 
estimates and also backward smoothness of trajectories. 
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Guide through the Proofs. 

We are not in a position to present complete proofs of the results. These are 
given in the cited references. Here, our aim instead is to orient the reader what 
are the main points one needs to go over in order to prove these resuhs. Special 
emphasis is given to more subtle technical details where "some" tricks -most of 
which recently introduced - need to be applied. It is our hope that an attentive 
reader will be able to grasp the essence of the proof, the tools required so that 
when guided to more specific reference he will be in a position to reconstruct a 
complete proof. 

Theorem 5.2 - Compactness 

Step 1: The corresponding system is gradient system. This follows 
from the fact that Lyapunov function V{t) = £{t) is bounded from above and 
below by the topology of the phase space. The latter results from dissipativity 
condition imposed on /(s) in Assumption 5.1 (see (3.4)). Moreover Vt{t) = 
implies {{g{wt),wt)) — hence Wt{x,t) = in 17 reduces the dynamics to a 
stationary set defined by (5.3) and (5.4). The set Af of stationary points is 
bounded due, again, to dissipativity condition imposed on /. This part of the 
argument is standard. 

Step 2: Asymptotic smoothness. The main difliculty in carrying the proof 
of asymptotic smoothness is double criticality of both the source and the damp- 
ing. This is the reason why previous approaches (based on splitting or squeezing) 
did require additional assumptions. To overcome this difficulty we use Theorem 
4.3. In preparation for this theorem we are using the equation 

zu-Az + aixMwt)-g{ut)]^f{w)~f{u) in O x M+ (5.14) 

with z = on r, written for the difference of two solutions z = w — u, where 
both w and u are solutions to the original equations confined to a bounded set 
B £ H, and also 

• energy identity (multiply equation (5.14) by Zj); 

• equipartition of energy (multiply equation (5.14) by z). 

In addition we recall dissipation relations for each solution u and w: 

f iigiut),ut))) + iigiwt),wt))<CB, > 0, (5.15) 

and also 

\\Mt)f + IWrnln + \\Mt)f + lk(t)ll?,o <Cb, vt > 0. 

Critical step consists of obtaining "recovery" estimate for the energy of z in 
terms of the damping and the source. The important part is that the "re- 
covery" estimate does not involve any initial data. This step is achieved by 
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using equipartition of the energy. The result of which is the following inequality 
[CL07a]: 

\tE,{T)+ f E,{t)<j [\\ztf + D{zt) + \G{zt,z)\]dt + -^{z,T), 
^ Jo Jo 

for T > 4, where E,it) = + l|Vz(t)||2, 

D{zt) = {{g{wt)-giut),zt)), G{zt, z) = {{g{wt) ~ g{ut), z)), (5.17) 

and 

■^{z.T)^- f {{f{w)-f{u),zt)) 
Jo 

+ f ds f dT{{f{w)-f{u),Zt))+[ \{{f{w)~f{u),z))\dt 
Jo Js Jo 

Our goal is to bound the right hand side of the above relation by terms that 
are (i) either bounded uniformly in T, or (ii) compact or (iii) small multiples of 
energy integral. 

From the energy relation for z we have that 

D{zt)dT <Cb+ f {{f{w) - f{u), Zt))dT. (5.18) 
Jo 

One can also wee from Proposition B.1.2 in [CLIO] that 
\\ztf <V + C^D{zt), V?7>0. 

Thus the terms Hztjp and D{zt) produce in (5.16) an admissible contribution. 

The last integral in the definition of ^ is compact. Thus the noncompact 
terms are the second one involving G and the first two integrals in definition of 

Since these two terms are similar, it suffices to analyze the main non-compact 
terms in (5.16) which are: 



Jo 



\G{zt,z)\dt and 



{{f{w)~ f{u),Zt))dT. 



These terms reflect double criticality of the damping-source exponents. We 
analyze these next. 

For the damping term G one uses split D, into Vti = {|u((x,t)| > 1} and 
the complement VL2 — Q,\VLi. With the use of critical growth condition for g 

and the fact that ||z|| < Cb wc obtain 



\g{ut)\\z\dx < 



\9{ut)f/'dx 



5/6 



< 



g{ut)utdx 



5/6 



<Gb I g{ut)utdx. 
Jq. 
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On 0,2 the damping is bounded -contributing to lower order terms. Thus we 
obtain for G{zt, z) 

G{zt,z)dt<CB / [igiut),ut) + {giwt),wt)]dt + CBT sup \\z{t)\\ 
Jo te[o,T] 

<Cb + CbT sup \\z{t)\\, 
te[o,T] 

whieh gives the inequahty in terms of the universal constant (independent on 
time T) and a lower order term ||z(t)||. 

For the source, the key point is the following decomposition that allows to 
prove sequential convergence of the functional ^ in Theorem 4.3: 

^iifiw)-f{u),zt)) 

[F{w{T)) + F{u{T)) - F{w{0)) - F{u{Qj)] 

^ m{u),wt)) + {{f{w),ut))]dt, 

where F{s) denotes the antiderivative of /. In line with Theorem 4.3 we con- 
sider these terms on sequences u = Wn,w = Wm converging weakly in (17) to 
a given element. We want to conclude that the corresponding functional con- 
verges sequentially to zero. Convergence on the first four terms is strong due 
to compactness while convergence on the last term is sequential with Wn ^ w 
weakly in H^[Vl) and Wmt ~^ Wt weakly in L2{0). Applying (5.16) gives 

E,{T)<- + ^ + C, [^{z, T) + LOT(z, T)] < e + C, [*(z, T) + LOT{z, T)] , 

for T = Ts large enough, where LOT{z, T) comprise of all compact terms. Now 
we use the fact that 

*(w„,Wm,T) = - / {{fiWn) - f(Wm),Wnt- Wmt)) 

ds / dT{{f{Wn) - f{Wm)-, W„t - Wmt)) 
J s 

is weakly sequentially compact - see the details in [CL07a, CL08a]. 

Step 3: completion. We apply Theorem 4.9 to deduce the results stated in 
Theorem 5.2 and Theorem 5.3. 

Theorem 5.3 -Structure -follows from asymptotic smoothness and gradient 
structure via Theorem 4.9. 

Theorems 5.4 and 5.7 - Finite Dimensionality and Regularity. The 

key tool is the quasi-stability inequality formulated in Definition 4.13. This 
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inequality, as in the case of asymptotic smoothness, should hold for the difference 
of two solutions z = w — u taken from a bounded set-say an absorbing ball B. 
As before, we apply two multipliers z and Zf. 

For the damping term G(zt, z) given in (5.17) we have the following estimate 
(see Proposition 5.3 in [CL07a] for details): 

Lemma 5.10 (Damping) For every e > there exists such that 



\G{zt,z)\dt<esnv\\z\\l,JT+ / [{{giwt),wt)) + iigiut),ut))]dT \ 

[0,T] \ Jo J 

+ C, [ {{g{ut)~g{wt),zt))dt. 
Jo 

The following decomposition of the source is also essential (see Proposition 5.4 
in [CLOTa]). 

Lemma 5.11 (Source) The following estimate holds true: 

i-T 



iif{u)^fiw),zt))<CB,TM' + ^ \\4i.n + Ce^B / Mt)\\i.nKmt, 

Jo Jo 

where K{t) = \\wt{t)\\'^ + ||ut(t)|p which under the condition g'{s) > m > 0, 
Vs, by (5.15) satisfies the following dissipativity property 

poo poo 

/ K{s)ds= / [\\wt{s)\\^ + \\ut{s)\\^]ds<CB. (5.19) 
Jo Jo 

By using dissipation properties in (5.15) and in (5.19) and combining with the 
results of Lemma 5.10 and Lemma 5.11 one obtains (see [CLOTa] for details) the 
inequality 

E^) <Ci [e-^'EM + MlMHl^in)] cxp £ X(s)ds| 
<C [e--'EM + Ml^^o,fM^^n)_ 

where 

Kit) = Kit) + iigiutit)),utit))) + iigiwtit)),wtit))). 
This leads to the desired quasi-stability inequality. 

Theorem 5.7 -Strong attractors - proof explores smoothness of backward 
trajectories on the attractor, see [KhalO]. 

Theorem 5.8-Decay rates. The details given in [CLOTa, Section 8]. 

Theorem 5.9 -Boundary Dissipation. The proof is more technical here. 
Let us point out the main differences. 

Compactness: 
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Gradient structure: In the boundary case the question of gradient structure 
has one additional aspect with respect to the interior damping. It involves a 
familiar in control theory question of unique continuation through the boundary. 
This property is known to hold due to various extensions of Holmgren-type 
theorems. This aspect of the problem is central in inverse problems when one 
is to reconstruct the source from the boundary measurements. Under suitable 
geometric conditions unique continuation results hold for the wave equation 
with a potential and solutions [Isa06, Rui92, EIN02]. 

Reconstruction of the energy: Since the damping is localized near the 
boundary, reconstruction of the energy involves propagation of the damping 
from a boundary into the interior. While in the case of interior damping equipar- 
tition of energy was sufficient for recovery of potential energy, this is not enough 
in the present case. Additional multiplier is needed which is referred as "flow 
multiplier" given by /iVw, where ft. is a suitable vector field. This multiplier -of 
critical order for the wave equation- is successful in reconstructing the energy 
of a single solution -provided time T is sufRciently large and the support of the 
damping on the boundary sufficiently large. However, when dealing with the 
differences of two trajectories (functions z), the method is no longer successful. 
There is a competition of two critical terms. A major detour is needed when 
one estimates the difference of two solutions. In the boundary case we have two 
critical multipliers /iVz and zt both of energy level: 

((/(") -/H,/^Vz)) and {{f{iL)- f{w),zt)). 

While this is not a problem with subcritical sources, it becomes a major issue 
in the case of critical sources. The remedy for the first term to this issue is by 
invoking Carleman's estimates with large parameter which allows for balancing 
the terms. This step is very technical-details can be found in [CLT09] . 

DiSSiPATiviTY INTEGRALS: Note that in this case (5.18) requires additional 
integration over the boundary in the right hand side. However, flux multiplier 
- again - with Carleman's estimates allows to get the estimate on this term. At 
the end of the day, with the use of restrictive growth conditions assumed on 
g (linear bound at infinity) one obtain reconstruction of the form as in (5.16 
without ll^tlp term. The estimate of G{z, zt) is reasonably straight due to linear 
bound hypotheses assumed. 

Functional ^E": The estimate for the compactness part of 'J is the same. 
It depends on the source only. This allows for the applicability of the same 
Theorem 4.3 to conclude compactness of attractors with boundary dissipation. 

Quasi-stability Consequences: Finite dimensionality and Smoothness 
of Attractors. Here the presence of the boundary damping presents another 
major challenge. This is at the level of the inequality in Lemma 5.11. The result 
stated there depends critically on finiteness of K{t). This is no longer valid in 
the pure boundary case, where dissipativity integral involves boundary integrals 
only. In that case one resorts to the analysis of " backward" smoothness on the 
attractor. This procedure consists of the following four steps. 
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Step 1: Prove quasi-stability inequality for solutions on the attractor near 
stationary points. This leads to consideration of negative time t — >■ — c». Such 
estimate is possible due to the fact that the system is gradient system and the 
velocities Wt{t) and ut{t) in (5.19) are "small" near equilibria. 

Step 2: Quasi-stability inequality provides additional smoothness of the trajec- 
tories at some negative time Tf < 0. This smoothness is propagated forward on 
the strength of regularity theorem. The consequence of this is that the elements 
in the attractor A are contained in H^(Jl) x H^{n). 

Step 3: The ultimate smoothness of attractors is achieved by claiming that the 
attractor is bounded in the topology of H^{il) x _ff^(f2). The above conclusion 
is obtained by exploiting compactness of attractors and finite net coverage. 

Step 4: Having obtained smoothness of attractor, one proceed to prove finite- 
dimensionality. This becomes a straightforward consequence of quasi-stability 
property on the attractor resulting from the aforementioned smoothness. 

It should be noticed that the method described above is fundamentally based 
on gradient property of the dynamics. The full details are given in [CLT09]. 

Remark 5.3 In the case of boundary damping (linearly bounded) there is an- 
other method that leads to compactness of attractors also in the presence of 
critical sources. This method relies on a suitable split of the dynamics [CEL02]. 
In such case there is no need for Carleman's estimates. However, properties 
such as regularity and finite dimensionality of attractors intrinsically depend on 
quasi-stability estimate. The latter requires consideration of the differences of 
two solutions. It is this step that forces Carleman's estimates to operate. In 
the subcritical case this is not necessary and more direct estimates provide full 
spectrum of results also in the boundary case, see [CEL04]. 

5.2 Von Karman plate dynamics 
Internal damping 

We consider system (2.4) with clamped boundary conditions (2.5) and discuss 
long time behavior under the following hypothesis. 

Assumption 5.12 We assume that Assumption 3.5 holds and, in addition, 
Qi G and g[{s) > m for \s\ > sq in the case a > and in the case a ~ 0: 
g G and g'{s) > m for \s\ > sq. 

We also want to be more specific about the nature of the source P{w). For 
sake of simplicity, we do not consider non-conservative force cases which lead 
to systems that are no longer of gradient structure. Instead we limit ourselves 
to conservatively forced models. We suppose that Assumption 3.6 holds with 
Pi(w) = 0, i.e., P(w) is a Frechet derivative of the functional — Ho(u') and the 
following property holds: there exist a e M and e > such that 

Ho(«;)-f a||u;||te,o>0, {{Po{w),w)) + a\\w\\l_,^,,>Q, Ww e H^in). (5.20) 
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Remark 5.4 In the specific example given in Remark 3.2 the functional Hq 
takes the form 

Tioiu) = -^{{[Fo,u],u))-{{p,u)), Fo e H^+\n) n H^n), peL^m. 

Thus, the assumption (5.20) is satisfied (in the case a > we can take less 
regular Fq). 

Under these conditions concerning P the energy function can be written as 
E{u, u,) = 1 [\\u,f + a\\Vu,f] + \\\I\Hn)f + I^o{u). 

Under the general Assumption 3.5 we have shown that von Karman system 
generates a continuous semiflow in the space Ha = -ffo(^) ^ ^a(i^), where 
Ha{Vl) is given by (3.18), i.e., Ha{^) ~ HQ{il) in the case a > and Ha{^) = 
L2(f^) in the case a = 0. This leads to a well-defined dynamical system {Ha, St). 

The existence of a global attractor under the assumptions specified for model 
(2.4) can be derived now from Theorem 4.3. The main long-time result reads 
as follows: 

Theorem 5.13 (Compactness) Let Assumption 5.12 he in force. Assume 
that P{w) = [Fo,w] +p, where Fq € H^+^{n) f] H^i^l) and p G L2in). Then 
the dynamical system {Ha, St) generated by equations (2.4) with the clamped 
boundary conditions (2.5) possesses a global compact attractor A in the space 
Ha = -ffo(fi) X 'Ha{^), where Ha is given by (3.18). 

In the case a > the proof of Theorem 5.13 follows from general abstract 
results for second order equations with subcritical source (see [CL08a]). In the 
case a ~ the force is critical and the proof relies on Theorem 4.3 in the same 
manner as in the case of wave equations. We refer to [CLIO] for details. 

Remark 5.5 One can easily obtain the same result for a more general class of 
nonlinear sources P{w) which can be axiomatized by (5.20) with an additional 
requirement that the map is compact. 

It is easy to verify (internal damping) that the energy £{u,ut) is a strict 
Lyapunov function for {Ha, St). Hence {Ha, St) is a gradient system (see Defi- 
nition 4.8 in Section 4.3), and the application of Theorem 4.9 yields the following 
result. 

Theorem 5.14 (Structure) Let Assumption 5.12 he in force. Assume that P 
in (2.4) satisfies (5.20). Let {Ha, St) be the dynamical system generated by equa- 
tions (2.4) with clamped boundary conditions Assume that {Ha, St) possesses a 
compact global attractor A. Let Af* be a set of stationary solutions to (2.4) (see 
Proposition 3.11). Then 

• A = M^{M), where M = {{w, 0) : w e AC} is the set of stationary point 
of {Ha, St) and M^{M) is the the unstable manifold Ai^{M) emanating 
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from Af which is defined as a set of all U € Ha such that there exists a 
full trajectory 7 = {U(t) = {u{t)]Ut{t)) : t € M} with the properties 

C/(0) = C/ and lim dist^^ (C/(t),7V) = 0. 

• The global attractor A consists of full trajectories 7 = {{u(t)\Ut{t)) : t G 
Ml such that 



t^±co y " weNt 

• Any generalized solution u{t) to problem (2.4) stabilizes to the set of sta- 
tionary points; that is, 

lim \\\ut{t)\\l^^+ M \\u{t)-w\\l\^0. (5.21) 

This theorem along with generic-type results on the finiteness of the number 
of solutions to the stationary problem allow us to obtain the following result. 

Corollary 5.15 Under the hypotheses of Theorem 5.14 there exists an open 
dense set TZq in L2{^) such that for every load function p € TZq the set Af of 
stationary points for (Ha, St) is finite. In this case A = ujIiA^"(zi), where 
Zi = (u'i;0) and Wi is a solution to the stationary problem Moreover, 

• The global attractor A consists of full trajectories 7 ~ {{u{t);ut{t)) : t S 
M} connecting pairs of stationary points; that is, any W € A belongs to 
some full trajectory 7 and for any "f G A there exists a pair {w-,w+} C 
Af* such that 

hm {\\ut{t)\\n^+\Ht)-w.\\l}^0 

and 

hm {\\utit)\\l^^ + \\uit)-w+\\l}^0. 

• For any (uq; ui) G Ha there exists a stationary solution w G such that 

hm {\\ut{t)\\l^ + \\u{t)-w\\l]^Q, (5.22) 

where u{t) is a generalized solution to problem (2.4) with initial data 
{uo]Ui) and with the clamped boundary conditions. 

Remark 5.6 We note that the generic class TZq of loads in the statement of 
Corollary 5.15 can be rather thin. The standard example (see, e.g., [Chu99, 
Section 2.5]) is the set TZq in the interval [0, 1] of the form 

7^o = {x e (0, 1) : 3rkeQ, |x - r^l < e2-^-i} , 
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where Q is the sequence of ah rational numbers. This set is open and dense in 
[0,1], but its Lebesgue measure less than e. However using analyticity of the 
von Karman force terms and the Lojasiewicz-Simon inequality we can obtain a 
result (see [Chullb]) on convergence for all loads p under some additional non- 
degeneracy at the origin type conditions concerning damping terms. By assum- 
ing this additional non-degeneracy condition (always satisfied when g'(0) > 
or in the hyperbolic case when g{s)s > ms^, \s\ < 1) the result of Corollary 
5.15 holds for all loads p € L2(f^)- For general discussion of Lojasiewicz-Simon 
method we refer to [HJ99, HJ09] and to the references therein. 

It follows from Assumption 5.12 (sec [CLIO], for instance) that there exists 
a strictly increasing continuous concave function Ho{s) such that 



in the case a > and < Ho{sg{s)) for s e R in the case a = 0. 

The above facts are used in formulation of the result on the rate of con- 
vergence to an equilibrium. 

Theorem 5.16 (Rates of convergence to equilibria) Assume that P sat- 
isfies (5.20) in (2.4) and the hypotheses above that guarantee the existence of 
a compact global attractor A for the dynamical system {Ha, St) generated by 
problem (2.4) with the clamped boundary conditions hold. 

If the set of stationary solutions consists of finitely many isolated equilibrium 
points that are hyperbolic,^ then for any initial condition y G Hq. there exists an 
equilibrium point e = (w; 0), w G Hq{Q,), such that 



where C and T are positive constants depending on y and e, [a] denotes the 
integer part of a and a{t) satisfies the following ODE, 



Q{s) = s - {I + Go) '(s) with Go{s)=ci{I + Ha) ' (css), (5.25) 



where positive numbers ci and C2 depend on y and e. If in Assumption 5.12 
So = 0, then the rate of convergence in (5.23) is exponential. 

Remark 5.7 It is also possible to state a result on the convergence rate with- 
out assuming hyperbolicity of equilibrium points (see [Chullb]). However as 
in the case described in Remark 5.6 this requires additional non-degeneracy 
type properties of the damping terms. Moreover, the lack of hyperbolicity is 
compromised by slowing down the decay rates. 

3ln the sense that the Hnearization of the stationary problem with the corresponding 
boundary conditions around every stationary solution has a trivial solution only. 



.s^ + sl < Ho{sig^{si) + S2.9j(s2)) for all si, S2 e R. 



\\Sty ~ e\\%^ < C ■ a {[tT-']) , t > 0, 



(5.23) 




(5.24) 
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The following result describes finite dimensionality of attractors. 

Tlieorem 5.17 (Dimension) Assume the hypotheses which guarantee the ex- 
istence of a compact global attractor A for the system (Ha, St) generated by 
problem (2.4) with the boundary conditions (2.5) hold (see Theorem 5.13) and 
a > 0. In addition, we assume that 

(i) In the case a ~ Q there exist to, Mq > such that 

< m < g' [s) < Mo[l + sg{s% s G R, (5.26) 

(ii) When a > Q we assume that (i) g{s) is either of polynomial growth at 
infinity or else satisfies (5.26), and (ii) there exists < 7 < 1 such that 
the functions gi satisfy the inequality 

< m < g[{s) < M[l + sg,{s)]'^, s G M, i 1, 2, (5.27) 
Then the fractal dimension of the attractor A is finite. 

Remark 5.8 One can see that in the case a > the condition imposed on gi 
in (5.27) arc valid if in addition to Assumption 5.12 we assume that that there 
exists TO. Ml, M2 > such that 

< g'ds) < Mi[l + \s\P~'], sg,is)>m\s\^P-^^/^ -M2, s&R, 

for some p> 1 (sec Remark 9.2.7 [CLIO]). 

The following assertion on regularity of elements in the attractor is a straight- 
forward consequence of Theorem 4.17. 

Theorem 5.18 (Regularity) Let the hypotheses of Theorem 5.17 be valid. 
Then any full trajectory 7 = {{u{t)]Ut{t)) : t € M} from the attractor A of 
the dynamical system (Ha, St) generated by problem (2.18) with the boundary 
conditions (2.5) possesses the properties 

u{t) e Cr{R;W), {ut{t),utt{t)) eCr{R;Ha), (5.28) 

where Cr means right- continuous functions and 

• in the case a>0: Ha = H^{^) x H^i^l) and W = H^{n) n H^i^). 

• in the case a = 0; Ha = H^{^) x ^2(1^) and W = H'^ifl) n H^{n). 
Moreover, 

• in the case a > 0: A C H^{Q) x H^{Q) and 

sup {\\uit)\\l + \\utit)\\l + \\uttit)\\l) <Ca<^, 
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• in the case a = 0: Ac H^{U) x H'^ 



(fl) and 



sup {\\u{t)\\l + \\ut{t)\\l + \\utt{t)\\l) <Ca<^. 



The validity of quasi-stability property for the system allows us to deduce, 
without an additional effort, another important property of dynamical system. 
This is existence of strong attractors that is, attractors in a strong topol- 
ogy determined by the generators of the (linearized) dynamical system. The 
corresponding result is formulated below. 

Theorem 5.19 (Strong attractors) We assume that the hypotheses of The- 
orem 5.18 are in force. 

Case a > 0: Let the rotational damping be linear; that is, gi{s) = gi ■ s for 
i = 1,2. Then the global attractor A is also strong: for any bounded set B from 
Hst = (H^i^l) n i?o(^)) X ^o(^^) we have that 



Case: a = 0: Let H^t = iH^{n) n H^i^l)) x H^{n). The global attractor 
A of the system (Ho^St) generated by the generalized solutions to problem (2.4) 
with clamped b.c. (2.5) is also strong, i.e. A is a global attractor for the system 
iHst,St). 

Moreover, in both cases a = and a > the global attractor A has a finite 
dimension as a compact set in Hgt- 

Proof. We apply the idea of Theorem 8.8.4[CL10]. The hypotheses of Theo- 
rem 5.18 guarantees that the system {Ha, St) is quasi-stable. I 

An important property of the dynamical system is an existence of exponential 
attractors, which attract trajectories at the exponential rate. The quasi-stability 
property allows us to deduce an existence of exponential attractors. The 
corresponding result is formulated below. 

Theorem 5.20 (Exponential attractors) Let {Ha, St) be the dynamical sys- 
tem generated by problem (2.4) and (2.5). 

Case a > 0: Under the condition (5.27) the dynamical system {Ha, St) 
has a (generalized) fractal exponential attractor A (see Definition 4.18j whose 
dimension is finite in the space H^{Q,) x W' , where W' is a completion of Hq{^1) 
with respect to the norm \\ • = ||(1 — aA) • j|_2- 

Case a = 0: Assume the validity of the hypotheses above which guarantee 
the existence of a global finite- dimensional attractor. Assume in addition that 



for some < 7 < 1. Then the system {HQ,St) has a (generalized) fractal 
exponential attractor A whose dimension is finite in the space (fi) x H^^{Q). 



lim sup dist/f ^ {Sty, A) = 0. 



\g{s)\ < C {1 + sg{s))'' , seR, 



(5.29) 
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Main ingredients for the proofs. The rotational case a > corresponds 
to the situation where von Karnian nonlinearity is compact, hence subcritical. 
This alleviates number of technical issues when proving validity of quasi- stability 
inequality which, in turn, leads to the existence of smooth and finite-dimensional 
attractors. Indeed, build in compactness turns critical integrals into lower order 
terms. 

We thus focus here on the more subtle non-rotational case, a = 0. Here, the 
main challenge is to establish quasi- stability inequality. Once this is proved, 
finite-dimensionality of attractors- Theorem 5.17, smoothness of attractor - 
Theorem 5.18 and existence of strong and exponential attractors -Theorem 5.19 
and also Corollary 5.20 follow from the abstract tools presented in Section 4.4. 

Thus, the core of the difficulty relates to the quasi-stability estimate obtained 
for a difference of two solutions z = u — w defined on Q = ^ x (0,r),S = 
r X (0, T). The solutions under consideration are confined to a bounded invariant 
set (attractor, for instance), so we can assume 

WMtW + hmin < R\ WMtW + wnmin < R'- 

To present the main idea we assume that P{u) = in (2.4). In this case the 
equation for z can be written as: 

d 

zu + A^z = f in n, z = 0, —z = on T, (5.30) 

on 

where / = ~a{g{ut) — g{wt)) + [v{u),u] — [v{w),w] . We have the following 
energy equality (satisfied for strong solutions) 

E4T) + dJ{z) = E4t)+ {n{z),zt)ndt, (5.31) 

where 

Mt) [V^t? + \Az\^]dx, n{z) = [v{u),u] - [v{w),w], 

and 




[g{ut) - g{wt)]ztdxdT. 



Our goal is to prove the following recovery inequality: 

E,{T) + f E,{t)dt < CrDI{z) + LOT{z), (5.32) 

"'0 

where 

LOT{z)<Cr sup |lz(i)||t.,a. 
te[o,T] 

Inequality (5.32), via reiteration on the intervals [mT, (m-t- l)r], leads to quasi- 
stability inequality. 
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To prove (5.32) one applies standard multipliers zt and z. After some cal- 
culations (as for the wave equation) one obtains 

E4T)+ [ E,{t)dt 
Jo 

< CrDI{z) + CrLOTIz) + Ct f iniz), zt)dt. (5.33) 

Jo 

The difficulty, however, is in handling the critical term 

/ {TZ{z),zt)ndt. 
Jo 

In order to obtain quasi-stability inequality, this critical term should be rep- 
resented in terms of the damping Dq{z), small multiples of the integrals of 
the energy, i.e., e Ez{t)dt and lower order terms (subcritical quantities) in a 
quadratic form. The key to this argument is "compensated compactness struc- 
ture" of the von Karman bracket which leads to the following estimate: 

I / {n{z),zt)dT\<CR max ||z(0||t. + C^fl / {\Wt\\ + \\wt\\)\\z\\ludt (5.34) 
Jo te[o,T] Jq 

The critical role is played by the presence in (5.34) of the velocities ||Mt||, ||wt|| 
which represent the damping and obey the estimate 

/•+00 

/ i\\ut\\'' + \\wt\\^)dt<CR (5.35) 
Jo 

Indeed, once (5.34) is proved and inserted into (5.33), then (5.35) along with 
Cronwall's inequality leads to the desired quasi-stability estimate. Thus, the 
key is to be able to prove (5.34). To this end the following "compensated 
compactness" decomposition of the von Karman bracket is used: 

in{z),z,)^~Qiz) + ^Pit) 

where 

Q{z) = {v{u) +v{w),[z,z]) -\Av{u + w,z)\'^, 
P{z) = -(wt, [w, u(z,z)]) - (wt, [w, u(z,z)]) 
~{ut + Wt, [z, v{u + w, z)]), 

where v{u,w) = A^^([u,w]). By using sharp Airy's regularity 

\v{u,v)\w2,^ < C||w||2,o||w||2,n 

given by Lemma 3.9, one obtains the estimates 

|QWI<C||z||t.,n and \P{z)\<CH{\\ut\\ + \\wt\\)\\z\\l^. 
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The above yields 



{niz),zt)dt < C{R)[Q{T)-Q{Q)]+C{R) I P{z)dt 

Jo 

< Cr f i\\ut\\ + \\wt\\)\\z\\l^ + CRina^ \\zml-e 
Jo *e[o.^l 

as desired for (5.34), hence for quasi-stability inequahty. 
Boundary damping 

We consider now the model in (2.18) with the boundary conditions in (2.19). 
This type of problems is often referred to as control problem with reduced 
number of controls. 

We impose the following hypotheses. 

Assumption 5.21 • Assumption 3.5(1,2) is in force, a G Looi^), a{x) > 
almost everywhere and P{w) = p G i2(^^)- 

• The boundary damping go G C^(M) is increasing, with the property go{0) = 
and there exist positive constants m, M , and sq such that 

< m < g'o{s) < M for \s\ > sq. (5.36) 

Case a > 0: 

We use the state space H = [-ff^(O) n H^i^l)] x H^{n). 

Using energy equality one can easily prove the following assertion. 

Theorem 5.22 Let Assumption 5.21 be in force and a > 0. Then 

• There exists i?* > such that the set 

WR^{y^{uo;ui)eH : 8{un,ui) < R} (5.37) 

is a nonempty invariant set with respect to semiflow St generated by equa- 
tions (2.18) with hinged b.c. (2.19) for every R > R^. Moreover, the set 
Wr is bounded for every R> R^ and any bounded set is contained in Wr 
for some R. 

• // in addition a{x) > almost everywhere in fl, the system {H,St) is 
gradient. 

• The set J\f of all stationary points of the semiflow St is bounded in H . 
Thus there exists i?*, > .such that M C Wr for all R > i?,,. 

We recall that 

J\f^{VeH : StV for all t > 0} . 
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and every stationary point W has the form W = (m; 0), where u = u{x) <E H'^iQ) 
is a weak (variational) solution to the problem 

A^u=[v{u),u]+p in rj; u = 0, Au = on T, (5.38) 

with function v{u) satisfying (2.6) with w ^ u, i.e., v{u) = A^^([u,u]). 

Theorem 5.23 (Compact attractors) Let a > and Assumption 5.21 be in 
force. Then 

• The restriction (W/j, S"*) o/ the dynamical system {H, St) on Wr given by 
(5.37) has a compact global attractor Ar C Wr for every i? > i?*, where 
i?* is the same as in Theorem 5.22. 

• //, in addition, a{x) > almost everywhere in Q, then there exists a 
compact global attractor A for the system {H,St). Moreover, we have 
A = A1"(A/'), where Ai'^{M) is the unstable manifold (see the definition 
in Section 4.3 ) emanating from the set M of equilibria for the semiflow 
St. 

• The two attractors Ar and A have finite fractal dimension provided rela- 
tion (5.36) holds for all s e M. Moreover, in this latter case the attrac- 
tors are bounded sets in the space H^{n) x H'^{n) and for any trajectory 
(u(t);ut(t)) we have the relation 

\Wtt{t)\\iM + \\ut{t)\\2M + \Ht)\\3,n < C, teR. 

An important ingredient of the argument is 

Proposition 5.24 (Observability estimate) Assume that a > and As- 
sumption 5.21 is in force. Let 

U{t)^{u{t)-ut{t))^Styi and W{t) ^ {w{t)-wt{t)) ^ Sty2 

be two solutions corresponding to initial conditions yi and y2 from the set Wr 
given by (5.37). Then there exist Tq > and constants Ci(T) and C2{R,T) 
such that 

2 

+ C2 [R, T) ■ LOT{z) (5.39) 
for any T >To, where z ~ u — w, 

E^t) - i [\zt\^ + a|Vzt|2 + |Az|2] dx, 
and the lower-order terms have the form 

LOT{z)^ sup \\z{t)\\1,,+ f \\zt{T)\\lndr. (5.40) 

0<T<T JQ 





d 


L 


on 
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We can obtain results on convergence rates of individual trajectories to 
equilibrium under the condition that the set Af of stationary points is discrete. 

As in the previous section we introduce a concave, strictly increasing, con- 
tinuous function Hq : M+ n- R_|_ which captures the behavior of go(s) at the 
origin possessing the properties 

iJo(O) = and < Ho{sgo{s)) for \s\ < 1. (5.41) 

We define a function Q{s) by relations (5.25) and consider the differential equa- 
tion 

^-fO(cr)=0, <>0, cr(0) = do e K+, (5.42) 
at 

which admits the global unique solution a(t) decaying asymptotically to zero as 
t — >■ oo. 

With these preparations we are ready to state our result. 

Theorem 5.25 (Rate of stabilization) Let the hypotheses of Theorem 5.23 
be valid with a{x) > a.e. Assume that there exist 7 > and sq > such 
that the interior damping g{s) satisfies the relation sg{s) > 75^ for \s\ < sq. 
In addition assume that problem (5.38) has a finite number of solutions. Then 
for any V ^ H there exists a stationary point E = (e; 0) such that StV —> E 
as t — > +00. Moreover, if the equilibrium E is hyperbolic in the sense that the 
linearization of (5.38) around each of its solutions has the trivial solution only, 
then there exist C, T > depending on V, E such that the following rates of 
stabilization 

\\StV -E\\h <Cai[tT-']), t>0, 

hold, where [a] denotes the integer part of a and a{t) satisfies (5.42) with ao 
depending on V, E Cz H (the constants Ci in the definition of Q also depend on 
V and E). In particular, if g'f^[Q) > 0, then 

\\StV -E\\H<Ce-''* 

for some positive constants C and w depending on V, E Cz H . 

Case a = and H = [H^ifl] n H^ifl)] x L-iin): 

We consider a model with nonlinear boundary dissipation acting via hinged 
boundary conditions which does not account a for regularizing effects of rota- 
tional inertia. Thus, the corresponding solutions are less regular than in the 
case of rotational models. 

The following analog of Theorem 5.22 is valid. 

Theorem 5.26 Let Assumption 5.21 hold and a = 0. Then 

• There exists i?* > such that the set 

Wr = = (uo; ui)eH : £{uo, ui) < R} (5.43) 

is a nonempty bounded set in H for all R> R^. Moreover, any bounded 
set B d H is contained in Wj^ for some R and the set Wr is invariant 
with respect to the semifiow St- 
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• // in addition a{x) > then the system (H, St) is gradient. 

Our main goal is to prove the global attractiveness property for the dynamical 
system {H,St). This property requires additional hypotheses imposed on the 
data of the problem. 

Assumption 5.27 The interior damping function g{s) is globally Lipschitz, 
i.e., \g{si) - g{s2)\ < M\si - S2I for all si, S2 £ K. 

This global Lipschitz requirement is due to the fact that only one boundary 
condition is used as a source of dissipation. 
Our main result is the following. 

Theorem 5.28 (Compact attractors) We suppose that Assumptions 5.21 
and 5.27 are in force. Then 

• For any R > there exists a global compact attractor Ar for the restric- 
tion {Wr, St) of the dynamical system {H, St) on Wr, where Wr is given 
by (5.43). 

• // we assume additionally that a{x) > a.e. in and g{s)s > for all 
s ^ 0, then there is Rq > such that Ar does not depend on R for all 
R > Rq. In this case A = Ar^, is a global attractor for {H,St) and A 
coincides with the unstable manifold M^[M) emanating from the set J\f 
of stationary points for St- Moreover, limt^+oo d-^stniStW^Af) = for 
any W € H. 

• The global attractors Ar^ and A are bounded sets in H^{fl) x H^(fl) and 
have a finite fractal dimension provided the relation in (5.36) holds for all 

s e R. 

From Theorem 5.28 we obtain the following corollary. 

Corollary 5.29 Let the hypotheses of Theorem 5.28 be in force. Assume that 
a{x) > a.e. in Q and g{s)s > for all s 7^ 0. Then the global attractor A 
consists of full trajectories j = {W{t) : t € M} such that 

lim distH{W{t),U) =0 and lim distH{W{t),JV) ^ 0. 

In particular, if we assume that equation (5.38) has a finite number of solutions, 
then the global attractor A consists of full trajectories 7 = {W^(i) : t € M.} 
connecting pairs of stationary points: any W € A belongs to some full trajectory 
7 and for any j d A there exists a pair {Z,Z*} C Af such that W{t) — >■ 
Z as t ^ —00 and W(t) — >■ Z* as t ^ +00. In the latter case for any V Cz H 
there exists a stationary point Z such that StV Z as t +00. 

In analogy with the previous cases one can also provide statements for the 
rate of stabilizations to equilibria. The statement (and the proofs) do not depend 
on the specific boundary conditions. However, in this case we do not need to 
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assume additional coercivity estimates for the interior damping. The point 
is that our basic Assumptions 5.21 and 5.27 are sufficient to conclude that 
generalized solutions are weak (variational). 

Main ingredients for the proofs: As in the case of interior damping, the 
rotational case a > is less involved. We shall thus focus on the non-rotational 
model a = where additional subtleties are present. Similar to the interior 
damping case, Theorem 4.3 provides the main tool for establishing asymptotic 
smoothness. For sake of avoiding repetitions, we shall mainly emphasize the 
parts of the proof that are more specific to the boundary damping. The proof is 
divided into several steps which are presented below. As before, we use notation 
g = X (0,T), E = r X (0,T) and assume that P{u) = in (2.18). 

We can assume that u and w are strong (smooth) solutions. Then the 
difference z = u — w solves the following problem 



where f = — 
are given by 



ztt ~ 
a{g{ut) 



Ah 



f 



.90 



in f2, 

[v{u),u] 

d 



-ut 



: 0, t/i on r, (5.44) 

[v{w), w] and the boundary conditions 

d 



90 



-Wt 



We have the following energy equality (satisfied for strong solutions) 

rT 



E,{T) + Di (z) = E,{t) + / {{TZ{z),zt))dt, 



(5.45) 



where 



and 



with 



J O 



72.(z) = [i;(u),w] — [v(w),w\ 



a[g{ut) - g{wt)]ztdxdT + DJ{z), 



go 



-Ut 



go 



-Wt 



-ZtdTdr. 



We also use the following assertion. 

Lemma 5.30 Let T > be given. Let (j) £ C^(]R) be a given function with 
support in [S,T — 5] , where 5 < T/A, such that < < 1 o,nd (j) = 1 on 
[2(5, T — 2(5]. Then any strong solution z to problem (5.44) satisfies the following 
inequality 



ES)mdt < Ci / ES)\4>'{t)\dt + C2 



fhV z(j)dxdt 



TZ{z)ztdxdt + C3 



d_ 

dn 

Bt[z) 



dE 



(5.46) 
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where the constants Ci do not depend on T and the boundary terms'^ Bt{z) is 
given by 



Bt{z)^ f 









2 


d d 


2 










it) ' 


+ 


dn dr 




<j)dj:+ [ 

Jo 


on 



(f)dt. (5.47) 



-i,r 



On the next step we eliminate the second- and third-order boundary traces 
on the boundary in the expression (5.47) for Bt{z). These second-order "su- 
percritical" boundary traces are due to the fact that the dissipation is allowed 
to affect the system via only one boundary condition. Traces of the second and 
third order (see (5.47)) are above the energy level, so these cannot enter the 
estimates. It is the microlocal analysis argument, again, that allows for the 
elimination of these terms. Handling of the boundary terms requires delicate 
trace estimates used already in the case a > 0. However the treatment of critical 
source (giving rise to the term TZ{z) ) is more troublesome now. The dissipa- 
tivity integral is supported on the boundary and not in the interior. Here are 
few details. For the trace result we use a more general trace estimate proved in 
[JL99] (see Proposition 1 and Lemma 4) valid for the linear KirchofF problem 
(note that the estimates in [JL99] are independent of the parameter represent- 
ing rotational forces). Thus, these estimates are applicable to both Kirchoff and 
Euler-Bcrnoulli (5.44) models. As a consequence we have the following 

Lemma 5.31 Let z be a solution to linear problem (5.44) with given f and ij) 
and Bt{z) be given by (5.47). Then there exist constants Ct > such that for 
any < rj < \ the following estimate holds 



Bt{z) < Ci.T I^IV-I' 



dn 



Zt 



dT. 



2.T 



ll^llc([0,T];_ff2-„(Q)) + ||2:t|li2([0,T];_ff-'i(n)) + / ll/ll-j;,n'^* 



The above calculations along with the Compactness Theorem 4.3 imply exis- 
tence of compact attractor -first statement in Theorem 5.28. As for smoothness 
and finite dimensionality, one needs to prove quasi-stability estimate -which 
amounts to the proof of the estimate in Proposition 5.24 -but with a = 0. This 
means that von Karman bracket is no longer subcritical and the term {TZ{z),Zt) 
needs to be estimated by appropriate damping on the boundary. To achieve 
this we proceed as in the internal case- by decomposing the bracket into Q{z) 
and P{z). However, in the boundary case the additional difficulty is due to 
the fact that dissipativity estimate (5.35) no longer holds. What one has, in- 
stead, is similar quantity on the boundary. In order to take advantage of this 
relaxed dissipation, one proceeds as in the case of wave equation by considering 



*They are not defined on the energy space. These are higher-order boundary traces of 
solutions. 
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(i) smoothing effect of backward trajectories on the attractor, (ii) propagating 
it forward and (iii) using the compactness of the attractor. This is technical 
part of the argument with full details given in [CL07b, Section 4.2]. 

Generalizations 

1. The case of nonconservative forces can also be considered. This is to say 
that the the force P{u) is not potential. In that case the system has no 
longer gradient structure and one needs first to establish existence of an 
absorbing set. This can be done under some additional restrictions on the 
damping, see [CLIO], where the case "potential operator + (7/),V)w" is 
considered. 

2. More general sources P{w) can be considered, including nonlinear sources. 
However, in general, one needs to make appropriate compactness hypothe- 
ses, see the Remark 5.5. 

3. Dissipation in other (e.g., free) boundary conditions, see [CLIO, CL04b, 
CL07b]. 

Open questions 

1. General theory of non-conservative forces which destroy gradient struc- 
ture. 

2. Localized dissipation, i.e., dissipation localized in a suitable layer near the 
boundary. 

3. Boundary damping problems without the light damping. This has to do 
with the existence of strict Lyapunov function and related unique con- 
tinuation across the boundary. While various unique continuation results 
from the boundary are available also for plates [AlbOO], the non-local na- 
ture of the von Karman bracket prevents Carleman's estimates [AlbOO] 
from applicability. 

5.3 Kirchhoff -Boussinesq plate 
Interior damping. Case a > 0: 

In this case our main results presented in the following theorem (for the proofs 
we refer to [CL08a, Chapter 7]). 

Theorem 5.32 Let a > and P{w) = A[w^] — q\w\^~^w for some p > and 
1>1. Assume that the damping functions g and G have the structure^ described 
in (2.8). Then problem (2.12) with a > and the clamped boundary conditions 
(2.5) generate a continuous semiflow St in the space H = H^{Q.) x H^{Vl). 



^see Remark 5.9(2) below. 
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Assume in addition that either 



{ 



m < 3, Q > 0, I < m, inf {a{x)} > 



or else 



{ 



TO < 3, Q = 0, inf {a{x)} > is sufficiently large 



} 



Then the semiflow St possesses a global compact attractor A. Moreover, if the 
constant Gi in (2.8) is positive, then 

• the fractal dimension of the attractor A is finite; 

• the system possesses a fractal exponential attractor (see Definition 4^.18) 
whose dimension is finite in the space Hq{VI)'><W , where W is a completion 
of Hq{Q) with respect to the norm \\ ■ \\w ~ II (1 ^ • 11-2- 

Remark 5.9 (1) As we see the additional potential term p|u|'^^u in the force 
P{w) allows us to dispense with a necessity of assuming large values for the 
damping parameter. 

(2) The requirements concerning the damping terms in Theorem 5.32 are 
not optimal. We choose them for the sake some transparency. The same results 
under much more general damping functions can be found in [CL08a, Chapter 7] . 

Interior damping. Case a ~ with linear damping function: 

We assume that g{s) = s and consider the source term of the form 



Then the dynamical system {H,St) generated by equations (2.12) with a = 
possesses a compact global attractor 2t. 

If a = 0, then the system {H,St) is gradient and thus 21 = A4"{Af) is 
unstable manifold in H emerging from the set Af of equilibria. 

We note that the system under consideration is not a gradient system when 
a 0. As a consequence, the proof of existence of global attractor can not be 
just reduced to the proof of asymptotic smoothness. One needs to establish first 
existence of an absorbing ball. For this we use Lyapunov type function of the 
form 



P{w) 



aA[w'^] ~ g\w\ 



w with (7, g > 0. 



(5.48) 




(5.49) 




(5.50) 
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with £ = J min{l, fc}, see [CL06b, CLll] for details in the case when a{x) = k 
is a constant. Here the energy functional has the form 



£(t)^h\w,m' + l\\^wm'' + \ I \vw{x,t)\Ux 

w{x,t)\Vw{x,t)\'^dx + / \w{x,t)\^+^dx. 

In order to complete the proof of the existence of global attractor, hence of 
Theorem 5.33, we use Ball's method (see Theorem 4.4). 

Let us take with e = fc/2 and denote ^'(i) = Vk/2{t). Using the fact 

that energy identity holds for all weak solutions one can easily see that the 
following equality is satisfied for this Lyapunov's function ^'(t): 



5t*(i) + fc«'W+ / {a{x) - k)\wtYdx = K{w{t)), 
Jn 

where 

K{w{t))^ i {a\Vw\^ - -{{a{x) - k)w)wtdx 
Jn 2 

-ka ( w\Vw\^dx-^ f \\/w\^dx - ^ / \w{x,t)\'+^dx. 
Jn 4 21 + 2 Jq 

It is clear that the term K{w{t)) is subcritical with respect to strong energy 
topology, therefore using the representation 

^{T)+ I e-'=(^-*)L(t)di = e~'='^5'(0) + / e-^'^'^-''^ K{w{t))dt, 



with 

L{t) = / {a{x)-k)\wt{t)\'dx, 
Jn 

we can apply Theorem 4.4 to prove the existence of global attractor. 

To establish the statement in the case ct = we note that in this case 
£"(!/;, w') is a strict Lyapunov function and thus the system is gradient. 

Now we consider long-time dynamics of strong solutions. 

Let Hst = n i?o(^)] X ^^o(^)- The argument given in the well- 

posedness section shows that the restriction Sf of the semiflow St on Hst is a 
continuous (nonlinear) semigroup of continuous mappings in Hst- 

Theorem 5.34 (Compact attractor for strong solutions) Let a = and 

a{x) = k > be a constant (for simplicity) . Then semiflow St^ has a compact 
global attractor in Hst- This attractor ^"^^ possesses the properties: 

• 21"* C 21; 
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21^* = A^"j(A^) is unstable manifold in Hst emerging from the set Af of 
equilibria; 



• has a finite fractal dimension as a compact set in H. 



St- 



Moreover, if we assume in addition I — 2m — 1 is an odd integer, then any 
trajectory (w{t);wt{t)) from the attractor possesses the property 

sup||w(")(OIU,a < C„,™ < oo, n,m = 0,1,2,..., (5.51) 

tgR 

where it;'")(t) ~ d"w{t). In particular, Qi'^* is a bounded set in C™(rj) x C™{fl) 
for each m = 0, 1, 2, . . and thus a''* C C°°{Ti) x C°°{Ti). 

The proof of Theorem 5.34 consists of several steps. The situation corre- 
sponding to gradient flows (we deal with the case ct = 0). 

Step 1 (dissipativity in Hst) by the "barrier" method: Due to dissipa- 
tivity in H and energy relation it is sufficient to consider dissipativity of strong 
solutions w{t) possessing properties 

/•oo 

\\wt{t)f + \\w{t)\\l<R,t>0, / \\wt{T)fdT<CB^<oo. (5.52) 

Jo 

Let w{t) be strong solution satisfying (5.52) and v{t) — Wt{t). Consider the 
functional 

G{v,v')^^{\\v'r + \\Awr + \\v\\' 

|VwnVw(a;)p +2|(Vw,Vw)h2|^ dxj, (5.53) 

and denote 

H{t) = H{v, vt) =ao + G{v, v^) + e^v, Vt), (5.54) 

where Oq and e are positive parameters. Then using Bresis-Gallouet-Sedenko 
type inequality (see Lemma 3.15 and Remark 3.4) and choosing e small enough, 
we obtain that 



dH{t) 
dt 



+ lH{t) < C},\\wtit)fH{t) In (1 + Hit)) + C 



2 



with positive 7 (see [CLll] for details). The above inequality leads to the 
dissipativity property by barrier's method. 

Step 2 (attractor in Hst)- To prove the existence of the attractor we need 
to establish asymptotic compactness of the system in C. For this we use again 
Ball's method. Let e = fc/2 and ao = in (5.54). Then one can see that 

^EMl + kH{t)^Kstiwit),Wt{t)), 
at 
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where 

Kstiw,v)^3 / (Vw,Vv)R2|Vvpda; 
Jn 

- I [f\w)v - v\vtdx ~ t: [ [/'(^)^' - v]vdx, 
Jn ^ Jn 

where /(s) = gi|s|'~^s. The term Kst{w,v) is obviously subcritical with respect 
to the topology in Hgt, therefore we can apply the same Ball's argument to 
prove the existence of global attractor 21^*. The inclusion SI'** C 21 is evident. 
One can also see that 21** = A^"t(7V). 

Step 3 (Finite dimension of the attractor in Hst)- To prove finite 
dimensionality of the strong attractor in Hst we use smoothness properties of 
trajectories and again the methods developed in [CL08a]. For this we need to 
establish the quasi-stability estimate in Hst- 

Proposition 5.35 (quasi-stability inequality in Hst) Let the hypotheses of 
Theorem 5.34 be in force. Let w{t) and w*{t) be strong solutions to the problem 
in question satisfying the estimates 

\\w{t)\U + \\wt{t)\\2 + \\wtt{t)\\ + \\w*{t)\U + \\wl{t)\\2 + \\wUt)\\ < R 

for all t CzM- and for some constant R > 0. Let w{t) ~ w{t) — w*{t). Then 

wmwl + WMml + Wmtm' < cie-^* (iiu;(o)ii^ + wwtml) 

+ C2 fe-'''^'--m\w{r)\\l + \\wt{r)r)dr 
Jo 

for all t > 0, where 7 > and Ci, C2 > may depend on R. 

The proof (see [CLll]) follows the strategy applied to von Karman evolution 
equations and presented in [CLIO, Sect. 9. 5. 3]. 

By Theorem 4.16 (see also Theorem 4.3 in [CLOSa]) the relation in Propo- 
sition 5.35 is sufficient to prove that the fractal dimension of 21^** as a compact 
set in Hst is finite. 

Boundary damping 

In the case of a > we can use the same methods (sec Section 10.3 in [CLIO] 
for the details) as for von Karman model (2.18) with boundary damping (2.19). 
The case of the boundary damping with a = is still open. 

6 Other models covered by the methods 
presented 

In this section we shortly describe several models whose long-time dynamics can 
be studied within the framework presented above. These include: (i) nonlocal 
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Kirchhoff wave model with strong and structural damping; (ii) finite dimensional 
and smooth attractors for thermoelastic plates; (iii) control to finite dimensional 
attractors in thermal-structure, flow-structure and fluid-plate interactions. We 
also consider dissipative wave models which arise in plasma physics. 

6.1 Kirchhoff wave models 

In a bounded smooth domain 57 C M'' we consider the following Kirchhoff wave 
model with a strong nonlinear damping: 

( -cr(||Vit||2)Aitt - (/)(||Vu||2)Au + /(u) /i(a;), x e Q., t > 0, 

{ (6-1) 

[ u|ao = 0, m(0) = uo, ut{0) = ui. 

Here A is the Laplace operator, a and are scalar functions specified later, 
f{u) is a given source term, h is a given function in L'^{Q,). 

This kind of wave models goes back to G. Kirchhoff (d = 1, (j>(s) = <j>Q + (j)is, 
(t(s) = 0, /(u) = 0) and has been studied by many authors under different sets 
of hypotheses (see, e.g., [Lio78] and also [Chul2] and the references therein). 
The model in (6.1) is characterized by the presence of three nonlinearities: the 
source, the damping and the stiffness. 

We assume that the source nonlinearity f{u) is a function possessing the 
following properties: /(O) = (without loss of generality), 

fif := liminf {5-^(5)} > -00, (6.2) 

and also (a) if d = 1, then / is arbitrary; (b) if = 2 then 

\f'{u)\ <C{1 + for some p > 1; 

(c) if d > 3 then either 

\f'{u)\<C{l + \u\P-^) with some 1 < p < = ^-i-^, (6.3) 

or else 

co\u\P'^ ~ ci < f {u) < C2 {I + \u\P-^) with some <p <p*,, (6.4) 

where p** = ( , Q > arc constants and S-|- = (s -I- |s|)/2. 

We note that the conditions above covers subcritical, critical and supercrit- 
ical cases. 

Under rather mild hypotheses concerning functions ip and a we can prove 
that the problem is well-poscd. This holds even without the requirement that (j) 
is non- negative (see the details in [Chul2]). However in order to study long-time 
dynamics of the problem (6.1) we need to assume that the functions a and (j) 
from C^(R+) are positive and either 

(f>{s)s — > +00 as s — > +00 and jjf = liminf {s^^ f{s)} > 0. 
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or else 

jiA, = liminf 4>{s) > and /irf,Ai + fJ.f > 0, 

where Ai is the first eigenvalue of the minus Laplace operator in 51 with Dirichlet 
boundary conditions. 

As a phase space we consider H = [i7g(fi)nLp+i(il)] x L2{^) endowed with 
partially strong^ topology: a sequence {{uq',Ui)} C "H is said to be partially 
strongly convergent to (uo;ui) S H if Uq — >■ uq strongly in i/Q(f2), Uq — >■ Mq 
weakly in Lp^i{Q) and — > ui strongly in L2{^) as n — > 00 (in the case when 
d < 2 we take 1 < p < 00 arbitrary). 

Under the conditions stated above (see [Cliul2]) problem (6.1) generates an 
evolution semigroup St in the space H. The action of the semigroup is given 
by the formula Sty = {u{t)]Ut(t)), where y = {uq;ui) G H and u{t) is a weak 
solution to (6.1). 

To describe the dynamical properties of St it is convenient to introduce the 
following notion. 

A bounded set 21 C H is said to be a global partially strong attractor for St 
if (i) 21 is closed with respect to the partially strong topology, (ii) 2t is strictly 
invariant (5*421 = 21 for all t > 0), and (iii) 21 uniformly attracts in the partially 
strong topology all other bounded sets: for any (partially strong) vicinity O of 
21 and for any bounded set B inH there exists = t^,{0, B) such that StB C O 
for all t > t^,. 

The main result in [Chul2] reads as follows: 

Theorem 6.1 Assume in addition that f'{s) > — c for all s Cz M. in the non- 
supercritical case (when (6.4) does not hold). Then the semigroup St given by 
(6.1) possesses a global partially strong attractor 21 in the space % . Moreover, 
21 C -Hi = [H'^{n) n Hl{n)] X Hl{n) and 

snx^ ( \\/^u(t)\\^ + \\yut{t)\\^ + \Wum\n+ \\Mr)fdT) < 

teR V Jt / 

for any full trajectory 7 = {(ii(t); Ut{t)) : t Cz M} from the attractor 21. We also 
have that 

2l = M+(7V), where Af={{u;0)'En:<j){\\A^/'^uf)Au + f{u)^h}. 

The attractor 21 has a finite fractal dimension as a compact set in the space 
[H^+''{n) n Hl{n)] X H'''{n) for every r < 1. 

To prove this result we first establish quasi-stability properties of St on different 
topological scales (see [Chul2] for details). 

In similar way we can also study a model with structural damping with a 
non-supercritical force of the form 

f utt+(^{\\A^'^uf)Au + a{\\A^l^uf)A'>ut+F(u)^Q, t > 0, 

{ (6-5) 



"It is obvious that the partially strong convergence becomes strong below supercritical 
level {Ul(Sl) C Lp+i{0)). 
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with 9 S [1/2, 1), where .4 is a hnear positive self-adjoint operator with domain 
'D[A) and with a compact resolvent in a separable infinite dimensional Hilbert 
space H . In this case we assume that the damping a and the stiffness factors 
arc positive functions and 

$(s) = / (t){C)d^ +00 as s -)■ +00. 
^0 

The nonlinear operator F is locally Lipschitz in the following sense 

\\A-'[F{u,)-F{u2)]\\ < L{e)\\A"\u,-U2)l yW^u,\\<Q, 

and potential, i.e., F{u) = 11' (u), where n(u) is a Ci-functional on V{A^/^), 
and ' stands for the Frechet derivative. We assume that n(u) is locally bounded 
on 2?(^^/^) and there exist 77 < 1/2 and C > such that 

77$(pi/2u||2) + U{u) + C>0, ueHi= V{A^I^) . 

This hypotheses cover critical case of the source F ^ but not supercritical. For 
details concerning the model in (6.5) we refer to [ChulO]. 

6.2 Plate models with structural damping 

Wc consider a class of plate models with the strong nonlinear damping the 
abstract form of which is the following Cauchy problem in a separable Hilbert 
space H: 

utt + D{u,ut) + Au + F{u) = 0, t > 0; u|t=o = "o, Mt|t=o "i- (6.6) 

In the case of plate models with hinged boundary conditions, A ~ (— A^)^, 
where A^j is the Laplace operator in a bounded smooth domain in with 
the Dirichlet boundary conditions. We have then that H = L2{^) and 

V{A) = {ue H-^{n) : M = Alt = on dil} . 

The damping operator D{u,ut) may have the form 

D{u,ut) = A [o-o(w)Aut] - div [ui{u, Vu)Vut] + g{u,ut), (6.7) 

where iTo(si), (Ti(si, S2, S3) and (?(si, S2) are locally Lipschitz functions of Si G K, 
i = 1,2, 3, such that ao{si) > 0, (Ji{si,S2, S3) > and g{si, 82)32 > 0. Also the 
functions cti and g satisfy some growth conditions. We note that every term 
in (6.7) represents a different type of damping mechanisms. The first one is 
the so-called viscoelastic Kelvin- Voight damping, the second one represents the 
structural damping and the term g{u,ut) is the dynamical friction (or viscous 
damping). We refer to [LTOO, Chapter 3] and to the references therein for a 
discussion of stability properties caused by each type of the damping terms in 
the case of linear systems. We should stress that the conditions concerning ai 
above allow to have a pure viscoelastic damping (i.e., we can take ai = and 
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5 = 0). However a similar results remains valid for the case when (Jq = and 
CTi is independent of Vu. In this case the presence (or absence) of the friction 
g(u,Ut) has also no importance for long-time dynamics. 

The nonlinear feedback (elastic) force F{u) may have one of the following 
forms (which represent different plate models): 

(a) Kirchhoff model: F{u) is the Ncmytskii operator 

u ~Ki ■ div{|Vu|'^Vu - /.ii|Vu|'"Vu} + K2 {Wfu - ^2|w|™w} ~p{x), 
where > 0, q > r > 0, Z > m > 0, /.i^ e M are parameters, p £ L2{i^). 

(b) Von Karman model: F{u) = -[u,T{u) + Fq] - p{x), where Fq G H'^{n) 
and p G L2{^) arc given functions, the von Karman bracket [u, v\ is given 
by (2.7) and the Airy stress function F{u) solves (2.6). 

(c) Bcrgcr Model: In this case the feedback force has the form 



Fiu) 



Au-p{x), (6.8) 



where k > and F e M arc parameters, p E L2{fl); for some details and 
references see, e.g., [Chu99, Chapter 4] and [CL08a, Chapter 7]. 

One can show that the system generated by (6.6) is quasi-stable. Thus methods 
presented in these notes apply. They allow to show the existence of global 
attractors of finite fractal dimension which in addition posses some smoothness 
properties. See [CKIO, CK12] for more details. 



6.3 Mindlin-Timoshenko plates and beams 

Let 17 C be a bounded domain with a sufficiently smooth boundary T. Lets 
v{x,t) ~ {vi{x,t),V2{x,t)) be a vector function and w{x,i) be a scalar function 
on Vl X ]R-|_ . The system of Mindlin-Timoshenko equations describes dynamics 
of a plate taking into account transverse shear effects (see, e.g., [LL88, Chap.l] 
and the references therein). This system has the form 

avtt + k ■ g{vt) ^ Av + k ■ {v + Vw) = -/o(w) + V^, [/i(w)] , (6.9) 

wtt + k ■ go{wt) - K ■ div{v + Vw) = -/2(w). (6.10) 
We supplement this problem with the Dirichlet boundary conditions 

vi{x, t) =V2ix, t) ^0, w{x,t) = on T x M+. (6.11) 

Here the functions Vi{x, t) and V2{x, t) are the angles of deflection of a filament 
(they are measures of transverse shear effects) and wix^t) is the bending com- 
ponent (transverse displacement). The vector fo(v) = (/oi(i'i, '^2); /o2(vi, i'2)) 
and scalar /i and /2 functions represents (nonlinear) feedback forces, whereas 
g{vi^V2) = (<?i(i^i)j <?2(i'2)) and go are monotone damping functions describing 
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resistance forces (with the intensity k > 0). The parameter a > describes 
rotational inertia of filaments. The factor k > is the so-called shear modulus 
(from mechanical point of view the limiting situation k — > +0 corresponds to 
plane strain and the case k — > +oo corresponds to absence of transverse shear). 
The operator A has the form 

Xi T~ 2 X2 2 XiX2 

l + f pa \-v a2 I q2 

2 '~'x\X2 2 a;i "T" 

where < z/ < 1 is the Poisson ratio. In ID case (dim 51=1) the corresponding 
problem models dynamics of beams under the Mindlin-Timoshenko hypothe- 
ses. For details concerning the Mindlin-Timoshenko hypotheses and governing 
equations see, e.g. [Lag89] and [LL88]. We also note that in the limit k — >■ -|-oo 
the system in (6.9) and (6.10) becomes Kirchhoff-Boussinesq type equation, see 
[Lag89, CL06a]. 

We refer to [CL08a] and [CL06a] for an analysis (based on quasi-stability 
technology) of long time behaviour of the Mindlin-Timoshenko plate under dif- 
ferent sets of assumptions concerning nonlinear feedback forces, damping func- 
tions and parameters. By the same method long-time dynamics in thermo- 
elastic Mindlin-Timoshenko model was studied in [Fas07, Fas09] . 

6.4 Thermal-structure interactions 

This problem has the form 

Utt - oiAe + A^u = B{u), u\t=o = uq, Ut\t=o = Wi, 

(6.12) 

Ot + TjAe + aAut = 0, 0\t=o^0o, 

where ^ is a linear positive self-adjoint operator with a compact inverse in a 
separable infinite dimensional Hilbert space H. The nonlinear term B{u) can 
model von Karman (as in (2.4)) and Berger (see (6.8)) nonlincarities. In this 
case the variable 9 represents the temperature of the plate. 

For application the methods presented we refer to [CL08b] and [CLIO, Chap- 
ter 11] in the von Karman case, see also [BC08] for the Berger nonlinearity. 

6.5 Structural acoustic interactions 

The mathematical model under consideration consists of a semilinear wave equa- 
tion defined on a bounded domain O, which is strongly coupled with the Berger 
or von Karman plate equation acting only on a part of the boundary of the do- 
main O. This kind of models, known as structural acoustic interactions, arise in 
the context of modeling gas pressure in an acoustic chamber which is surrounded 
by a combination of hard (rigid) and flexible walls. 

More precisely, let O C M'^ be a bounded domain with a sufficiently smooth 
boundary dO. We assume that dO ^QUS, where 17 n S" = 0, 

V,C {x= {xi;x2; 0) : x' = {xi;x2) G M^} 



A = 
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with the smooth contour T = dfl and S is a surface which hes in the subspace 
Kil = {x3 < 0}. The exterior normal on dO is denoted by n. The set f2 is 
referred to as the elastic wall, whose dynamics is described some plate equa- 
tion. The acoustic medium in the chamber O is described by a scmilinear wave 
equation. Thus, we consider the following (coupled) PDE system 

zu+g{zt)-Az + f{z)^0 in Ox(0,T), 

dz , . dz „ . , (6.13) 



^ -0 on S X (0,T), —^avt on r2x(0,T), 
on on 



and 



(6.14) 



j vtt + b{vt) + A'^v + B{v) + Pzt\n = in nx{0,T), 

\ v = Av = on dnx (0,T), 

endowed with initial data 

z{0,-) = z°. zti0r)^z^ in fl, v(0,-)=w°, Wt(0, •) = u^ in n. 

Here above, g{s) and b{s) are non-decreasing functions describing the dissipation 
effects in the model, while the term /(z) represents a nonlinear force acting 
on the wave component and B(v) is (nonlinear) von Karman or Berger force; 
a and /3 are positive constants; The part 5* of the boundary describes a rigid 
(hard) wall, while is a flexible wall where the coupling with the plate equation 
takes place. The boundary term /Sztjn describes back pressure exercised by the 
acoustic medium on the wall. 

Well-posedness issues of the abstract second order system, which is a par- 
ticular case of the one studied in [Las02, Sect. 2.6] (see also [CLIO, Chapter 
6]). 

Long-time dynamics from point of view of quasi-stable systems were inves- 
tigated in [CLIO, Chapter 12] in the von Karman case and in [BCL07] for the 
Berger nonlinearity. We also note that flow structure model (6.13) and (6.14) in 
combination with thermoelastic system (6.12) was studied in [BC08] and [CLIO, 
Chapter 12]. 



6.6 Fluid-structure interactions 

Our mathematical model is formulated as follows (for details, see [CRll]). 

Let O C M"^, C and the surface S be the same as in the case of 
the model in (6.13) and (6.14). We consider the following linear Navier-Stokes 
equations in O for the fluid velocity field v = v{x, t) = {v^{x, t); v'^{x, t); v^{x, t)) 
and for the pressure p{x,t): 

Vt~vAv + Vp = Gf and divv:=0 in Ox(0,+oo), (6.15) 

where i/ > is the dynamical viscosity and G/ is a volume force. We supplement 
(6.15) with the (non-slip) boundary conditions imposed on the velocity field 
V = t;(.T, t): 

v = onS'; w = (w^; w^; w^) = (0; 0; m*) on rj. (6.16) 
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Here u = u{x,t) is the transversal displacement of the plate occupying S7 and 
satisfying the following equation: 

utt + A'^u + F{u) = p\n in r2x(0,c»). (6.17) 

The nonlinear feedback (elastic) force F{u) as above may have one of the forms 
(Kirchhoff, Karman or Berger) which are present for plates models with struc- 
tural damping (see (6.6)). We also impose clamped boundary conditions on the 
plate 

Du 

u\an=jr =0 (6-18) 
and supply (6.15)-(6.18) with initial data of the form 

w(0) = wo, u{0)^uo, ut{0)=ui, (6.19) 
We note that (6.15) and (6.16) imply the following compatibility condition 

u{x',t)dx' = const for all t>0, (6.20) 

which can be interpreted as preservation of the volume of the fluid (we can 
choose this constant to be zero). 

It was shown in [CRll] that equations (6.15)-(6.20) generates an evolution 
operator St in the space 

y. = |(wo;uo;wi) 'E X X -ffo(fi) X L2(ri) : ivo,n) = = ui on fij , 

where 

X = {v^ {v^;v^;v^) £ [^2(0)]^ : divw = 0; (v^n) = on S} , 

and L2(f^) (resp. HQ{fl)) denotes the subspace in L2{^) (resp. mHQ{il)) con- 
sisting of functions with zero averages. 

This evolution operator St is quasi-stable and thus possesses a compact 
global attractor, see [CRll]. We emphasize that we do not assume any kind 
of mechanical damping in the plate component. Thus this results means that 
dissipation of the energy in the fluid due to viscosity is sufficient to stabilize the 
system. 

In a similar way (see [Chulla]) the model which deals only with longitudinal 
deformations of the plate neglecting transversal deformations can be also con- 
sidered (in contrast with the model (6.15)-(6.19) which takes into account the 
transversal deformations only). This means that instead of (6.16) the following 
boundary conditions arc imposed on the velocity fluid field: 

w = on 5; v = {v^ ;v^;v^) = {ul;u^;0) on fi, 

where u = {u^(x,t)]u'^(x,t)) is the in-plane displacement vector of the plate 
which solves the wave equation of the form 

utt-Au-V[diYu] + iy{vl -vl )\^^^o + f{u) = in Q; = on T. 
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This kind of models arises in the study of blood flows in large arteries (see the 
references in [Gro08]). 

One can also analyze the corresponding model based on the full Karman shell 
model with rotational inertia (sec [CR12]). In this case to obtain well-posedness 
we need to apply Sedenko's method. The study of long-time dynamics is based 
on J. Ball's method (see Theorem 4.4). 

6.7 Quantum Zakharov system 

In a bounded domain f2 C M'', c? < 3, we consider the following system 
ntt - A (n + |£;p) + h'^A'^n + ant = fix), x e fl, t > 0, 

(6.21) 

iEt +AE- h'^A^E + i-fE - nE = g{x), x £ t > 0. 

Here E{x,t) is a complex function and n{x,t) is a real one, /i > 0, a > and 
7 > are parameters and /(x), g{x) are given (real and complex) functions. 

This system in dimension d = 1 was derived in [GHG05], by use of a quan- 
tum fluid approach, to model the nonlinear interaction between quantum Lang- 
muir waves and quantum ion-acoustic waves in an electron-ion dense quantum 
plasma. Later a vector 3D version of equations (6.21) was suggested in [IIS09]. 
In dimension d = 2, 3 the system in (6.21) is also known (see, e.g., [SSS09] and 
the references therein) as a simplified "scalar model" which is in a good agree- 
ment with the vector model derived in [IIS09] (see a discussion in [SSS09]). 

The Dirichlet initial boundary value problem for (6.21) is well-posed and 
generates a dynamical system in an appropriate phase space (see [Chullc]). 
Using quasi-stability methods one can prove the existence of a finite-dimensional 
global attractor, for details we refer to [Chullc]. 

We also note that in the case /i = we arrive to the classical Zakharov system 
(see [Zak72]). Global attractors in this case were studied in [Fla91, GM98] in the 
ID case and in [GS05] in the 2D case. In the latter case to obtain uniqueness the 
Sedenko method was applied and Ball's method was used to prove the existence 
of a global attractor. 

6.8 Schrodinger— Boussinesq equations 

The methods similar to described above can be also applied in study of qualita- 
tive behavior of the system consisting of Boussinesq and Schrodinger equations 
coupled in a smooth (2D) bounded domain C M^. The resulting system takes 
the form: 

wtt + liwt +A^w-A (fiw) + \E\'') = gi{x), (6.22a) 

iEt + AE ~wE + ij2E ^ g2{x), x e D,, t > 0, (6.22b) 

where E(x, t) and wlx, t) are unknown functions, E{x, t) is complex and w{x, t) 
is real. Here above 71 and 72 are nonnegative parameters and .91(2;) and g2{x) 
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are given (real and complex) L2-functions. We equip equations (6.22) with the 
boundary conditions 

w\dn = Aw\on = 0, E\an = 0, (6.23) 

and with the initial data 

wt{x,0)=wi{x),w{x,0)=wo{x), E{x,0) ^ Eo{x). (6.24) 

Long-time dynamics in this system was studied in [CSll] by the methods de- 
scribed above under the following hypotheses concerning the (nonlinear) func- 
tion f: f GC^ (R), /(O) = and 

3ci, C2 > : F{r) = f fi^d^ > -cir^ - C2, V|r| > ro, (6.25) 

3M>0,p>l: \f'{s)\<M{l + \s\P-'^), seR. (6.26) 

We also note that the ideas which were used in the study of the Kirchfoff- 
Boussinesq system (see (2.12) with a = can be also applied to the following 
Schrodinger-Boussinesq-Kirchfoff model: 

wtt - Awt + A^w - div {|Vu|^Vu} ~ A {w'^ + \E\'^) ^ gi, 

iEt + AE -wE + i-i2E = 52, a; e il, i > 0, 

equipped with the boundary conditions (6.23) and initial data (6.24). For details 
we refer to [CS12]. 
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